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ABSTRACT 



Aims. We review the present observational knowledge on the spatial distribution and the physical state of the different 
(molecular, atomic and ionized) components of the interstellar gas in the innermost 3 kpc of our Galaxy - a region 
which we refer to as the interstellar Galactic bulge, to distinguish it from its stellar counterpart. 

Methods. We try to interpret the observations in the framework of recent dynamical models of interstellar gas flows in 
the gravitational potential of a barred galaxy. 

Results. Relying on both the relevant observations and their theoretical interpretation, we propose a model for the 
space-averaged density of each component of the interstellar gas in the interstellar Galactic bulge. 
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1. Introduction 

The state of knowledge on the Galactic center (GC) en- 
vironment prevailing one decade ago was reviewed inde- 
pendently by Morris & Serabyn (1996) and by Mezger et 
al. (1996). To summarize the main points relevant to our 
work, the Galactic bulge (GB) is the region of our Galaxy 
inside Galactocentric radius r ~ 3 kpc - a radius which 
roughly corresponds to the inner boundary of the Galactic 
disk molecular ring. The outer parts of the GB, outside 
r ~ 1.5 kpc, are largely devoid of interstellar gas. Atomic 
gas is confined to a noticeably tilted layer extending (in 
projection) out to r 1.5 kpc (Burton & Liszt 1978; Liszt 
& Burton 1980). This Hi layer is sometimes referred to as 
the Hi nuclear disk (Morris & Serabyn 1996), but we find 
this denomination misleading, as we think that the term 
"nuclear" should be reserved to objects and regions with 
smaller Galactic radii (e.g., r < 300 pc, the approximate 
radius of the nuclear bulge). For want of any better term, 
we will refer to this Hi layer as the Hi GB disk. The GB 
disk also includes significant amounts of H2, whose gen- 
eral distribution and kinematics were argued to be similar 
to those of Hi (Liszt & Burton 1978; Burton & Liszt 1992; 
but see also Sanders et al. 1984; Combes 1991 for a different 
viewpoint). However, molecular gas tends to concentrate in 
the so-called central molecular zone (CMZ), an asymmetric 
layer of predominantly molecular gas extending (in projec- 
tion) out to r '--^ 200 pc - more exactly, r ~ 250 pc at posi- 
tive longitudes and r ~ 150 pc at negative longitudes. The 
CMZ itself contains a ring-like feature with mean radius 
~ 180 pc, now known as the 180-pc molecular ring, and, 
deeper inside, a thin sheet populated by dense molecular 
clouds, known as the GC disk population or the GC molec- 
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ular disk (Bally et al. 1987, 1988). The observed kinematics 
of the 180-pc molecular ring indicate strongly non-circular 
motions, which in turn suggest that either the ring is radi- 
ally expanding (Kaifu et al. 1972; Scoville 1972; hence the 
historically often used denomination of "expanding molec- 
ular ring") or, more likely, that the gas travels along highly 
elongated orbits (e.g.. Bally et al. 1988; Binney et al. 1991). 

A number of significant advances were made in the last 
decade, which contributed to improving our knowledge and 
understanding of the interstellar GB. Nonetheless, several 
aspects remain ambiguous or controversial, and important 
pieces of the puzzle are still missing, mainly due to the 
severe interstellar extinction along the line of sight to the 
GC, to the existence of strongly non-circular motions and 
to the lack of accurate distance information. Thus, to date, 
no complete and fully consistent picture has emerged from 
the vast body of existing observational data. 

In this paper, we put together a model for the spa- 
tial distribution of interstellar gas in the Galactic region 
r < 3 kpc, based on the observational results summed up 
by Morris & Serabyn (1996) and Mezger et al. (1996), on 
more recent observational studies of dust thermal emission, 
CO line emission. Hi 21-cm line emission and absorption 
and pulsar dispersion measures, and on recent theoretical 
investigations of interstellar gas dynamics near the GC. Our 
model complements that developed by Ferriere (1998, 2001) 
for the Galactic disk outside r ~ 3 kpc. In section [2l we re- 
view the observations pertaining to the spatial distribution 
and the physical state of interstellar gas in its molecular, 
atomic and ionized forms. In section [3l we discuss a few 
dynamical models of interstellar gas subject to the barred 
gravitational potential of the Galaxy, and we make the link 
between the theoretical predictions of these models and the 
observational facts reviewed in section [21 In section 21 we 
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reconcile as well as possible the disparate observational and 
theoretical results presented in the preceding sections, and 
we incorporate them into a new, consistent model for the 
space-averaged density of interstellar gas in the Galactic 
region r < 3 kpc. In section [5l we conclude our study. 

For consistency with Ferriere (1998, 2001) and in line 
with most papers discussed below, the Sun is assumed to 
lie at a distance = 8.5 kpc from the GCQ even though 
recent work favors a somewhat smaller value. All spatial dis- 
tributions will be given either as functions of Galactocentric 
cartesian coordinates {x,y,z), with the x-axis pointing to- 
ward the Sun, the y-axis in the direction I = +90° and the 
z-axis toward the North Galactic Pole (NGP), or as func- 
tions of Galactocentric cylindrical coordinates (r, 9, z), with 
9 increasing in the direction of Galactic rotation, i.e., clock- 
wise as seen from the NGP (see Figure[T]) . Toward the vicin- 
ity of the GC, y ~ (150 pc) (^1°) and z ~ (150 pc) (6/1°), 
where as usual I and b denote Galactic longitude and lat- 
itude, respectively. Finally, the projected (onto the plane 
of the sky) horizontal distance from the GC is given by 
r^ = \y\ 0,(150 pc) (|Z|/1°). 



2. Observational overview 

In this section, we successively review the current observa- 
tional status of the total interstellar matter (section l2.ip . 
the molecular gas (section [521), the atomic gas (section [53]) 
and the ionized gas fsection l^^ . For each, we provide quan- 
titative estimates for the space-averaged density and the 
mass of interstellar hydrogen in subregions of the GB. The 
density and mass of the total interstellar matter (including 
helium and metals) immediately follow, once the interstel- 
lar elemental abundances in the GB are known. 

Various lines of evidence, primarily from Hii regions 
(Shaver et al. 1983; AfHerbach et al. 1997), planetary 
nebulae (Maciel & Quireza 1999) and early B-type stars 
(RoUeston et al. 2000), converge to point to the existence 
of an inward metallicity gradient ~ (0.06 — 0.07) dex kpc~^ 
in the inner Galactic disk. If extrapolated in to the GC, 
this metallicity gradient implies an increase in metallicity 
by a factor ~ 3 — 4 between the Sun and the GC. On the 
other hand, direct abundance measurements near the GC 
yield mixed results : while some confirm the high metallicity 
predicted for the GC region (e.g., Smartt et al. 2001), oth- 
ers indicate instead a close-to-solar metallicity (e.g., Carr 
et al. 2000; Najarro et al. 2004). Here, we adopt the inter- 
mediate, and rather conservative, value Zqq = 2 for the 
GC metallicity (as in Launhardt et al. 2002, and consistent 
with Sodroski et al. 1995; see below). Furthermore, since 
the measured IIc+/II+ ratio does not exhibit any significant 
radial gradient (Shaver et al. 1983), we adopt Yqc = Yq 
for the GC helium fraction. Altogether, with a helium-to- 
hydrogen mass ratio of 0.4 and a metal-to-hydrogen mass 
ratio of 2 X 0.0265 (see, e.g., Anders & Grevesse 1989; 
Dappen 2000), the conversion factor from interstellar hy- 
drogen masses to total interstellar masses is 1.453. 



All our estimates can easily be rescaled to any other value 
of tq : distances scale as tq, surface densities as Tq, volume 
densities as Tq and masses as Tq. 




Fig. 1. The (x, y, z) and (r, 9, z) Galactocentric coordinate 
systems used in the present paper: (a) two-dimensional view 
looking down from the NGP; (b) full three-dimensional 
view. Note that both systems must be left-handed if one 
wants to stick to the usual conventions that z has the same 
sign as b (i.e., positive/negative in the northern/southern 
Galactic hemisphere) and that 9 increases in the direction 
of Galactic rotation (i.e., clockwise about the z-axis). 



2.1. Total interstellar matter in the innermost 500 pc 

The large-scale spatial distribution of interstellar matter in 
the CMZ was derived by Launhardt et al. (2002), based 
on IRAS and COBE/DIRBE far-infrared maps of the cen- 
tral kpc of the Galaxy. They converted the optically thin 
240 /im dust emission map into a hydrogen column den- 
sity map, using for the dust temperature an empirical color 
temperature proportional to the 60-to-140 /im flux ratio 
and for the H-to-dust mass ratio the standard value near 
the Sun, (Mn/Afd)© = 110, divided by the relative metal- 
licity, (Zgc/Zq) = 2. 

In contrast to previous authors, Launhardt et al. (2002) 
defined the CMZ as the entire region interior to r ^ 500 pc, 
and they divided it into the inner CMZ, physically associ- 
ated with the stellar nuclear bulge (r < 230 pc) and heated 
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Region of space Radial range Hydrogen mass Reference 



Inner CMZ : 


inner disk 


~ 0- 


120 


pc 




~ 4 X 10'' Mq 


Launhardt et al. (2002) 




outer torus 


~ 130 


pc- 


- 230 


pc 


~ 1.6 X lO'^ Mq 




Outer CMZ 


I > 0° 
I < 0° 


~ 230 




- 500 


pc 


~ 2.9 X lO'^ Mq 

~ 1.1 X 10'^ Mq 


!) 
!) 


Entire CMZ 




~ 0- 


500 


pc 




~ 6 X 10'' Mq 


!) 



by its stars, and the outer CMZ, containing dust too cold 
to be possibly heated by stars from the nuclear bulge. 

Extinction arguments as well as comparisons with high- 
resolution submm continuum maps led Launhardt et al. 
(2002) to conclude that interstellar matter in the nuclear 
bulge is extremely clumpy. According to them, > 90% of 
the interstellar matter would be trapped in small, compact 
molecular clouds occupying only a few % of the interstellar 
volume. The average hydrogen density inside these clouds 
would, therefore, be ~ a few 10'^ cm"'', with possibly much 
higher values in the cloud cores. The remaining < 10% 
of the interstellar matter would form a diffuse, homoge- 
neously distributed intercloud medium, Launhardt et al. 
(2002) further divided the inner CMZ into a warm inner 
disk with radius ~ 120 pc and a cold outer torus centered 
on the 180-pc molecular ring and extending radially be- 
tween ~ 130 pc and 230 pc. The whole disk/torus struc- 
ture has a FWHM thickness ~ 45 pc and a total hydrogen 
mass 2.0 X 10"^ Mq, with - 4 X 10^ Mq in the inner disk 
and ~ 1.6 x 10^ Mq in the outer torus. The space-averaged 
density of hydrogen nuclei can be approximated bj0 



('^H)disk('', z) = (72 cm ^) exp 



r — 70 pc 



X exp 



(1) 



in the inner disk and 

("-h) torus (?", ^) = (106 cm"^) exp 



r- 180 pc 

z 



X exp 



(2) 



in the outer torus, with L = (100 pc)/(2 {\n2y/*) ~ 55 pc 
and H = (45 pc)/(2 (ln2)i/i '^) ~ 29 pc. Note that strictly 
speaking. Equation [1] implies that the inner disk extends 
radially between ~ 20 pc and 120 pc, whereas the actual 
matter distribution is likely to continue in to the origin; 
but since the region r < 20 pc contributes very little to the 
total mass, this approximation is perfectly valid. 

Launhardt et al. (2002) also found large amounts of in- 
terstellar dust too cold to belong to the nuclear bulge, al- 
though its projected distribution onto the plane of the sky 
would place most of this dust inside the nuclear bulge. The 
associated interstellar hydrogen has a mass ~ 2.9 x 10^ Mq 
at positive longitudes (0?5 < / < 4°) and 1.1 x 10^ Mq 



^ The space-averaged densities given by Launhardt et al. 
(2002) must be divided by a factor of 3.5 in order to match the 
masses and column densities quoted in their paper (Launhardt, 
private communication) . 



at negative longitudes (—4° < I < —2°). It could be this 
interstellar matter, located in the outer CMZ, that would 
be responsible for the observed asymmetry of the CMZ. 
Altogether, the total mass of interstellar hydrogen enclosed 
within the central kpc of the Galaxy is '--^ 6 x 10^ M© (see 
Table HI). 

with an average hydrogen density 10 cm~^, and the 
strong UV radiation field from the numerous high-mass 
stars embedded in the nuclear bulge would cause this thin 
intercloud medium to be both warm and ionized. As will 
become apparent in the course of this section, the two above 
density estimates are compatible with more direct density 
estimates for molecular gas (see section and for ionized 
gas (see section \2A^ . 

2.2. Molecular gas 

Radio emission lines of interstellar molecules (excluding H2, 
which possesses no permitted lines at radio frequencies) 
constitue one of the best tools to probe interstellar molecu- 
lar gas in remote Galactic regions, and chief amongst them 
is the ^^CO (J = 1 ^ 0) line at 2.6 mm. Unfortunately, 
the ^^CO (J = 1 ^ 0) line as a tracer of molecular gas is 
plagued by the considerable uncertainty in the ^^CO-to-H2 
conversion factor, A"coj which relates the H2 column den- 
sity, Nb_2, to the velocity-integrated intensity of the ^^CO 
line, Wi2cQ, through Nn^ = Xqo Wi^qq. 

Sodroski et al. (1995) estimated the value of Xco near 
the GC by combining the Goddard-Columbia surveys of 
i^CO (J = 1 0) emission with COBE/DIRBE observa- 
tions at 140 /im and 240 /im. To deduce the H2 column 
density from the far-infrared data, they assumed that the 
gas-to-dust mass ratio is inversely proportional to metallic- 
ity, Z, and that Z is ~ 1.5—3 times higher near the GC than 
in the inner Galactic disk (2 kpc < r < 7 kpc). Accordingly, 
they found that Xco is ~ 3—10 times lower near the 
GC than in the inner Galactic disk, which, together with 
Xco 2.2 X 10^*^ cm~^ km~-^ s in the inner Galactic 
disk, leads to Xqo (2 — 7) x 10"^^ cm~^ K^^ km~^ s near 
the GC. The corresponding H2 mass inside the longitude 
range — 1?5 < I < 3?5 (projected half-size ~ 375 pc) is then 
- (2-6) X 10^ Mq. 

Arimoto et al. (1996) studied the radial dependence of 
the conversion factor in more detail, by comparing the CO 
luminosities to the virial masses of a large number of giant 
molecular clouds situated at various distances from the GC. 
In this manner, they obtained the relation (rescaled to r© = 
8.5 kpc) 



Xco{r) ~ (9 X lO'-" cm"^ K"' km"' s) cxp 



7.1 kpc 



(3) 

Note in passing that Equation [3] implies Xco {fo ) ~ 3 x 
10^" cm-2 km-i s. 
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Table 2. Estimated values of the ^^CO-to-H2 conversion factor, Xqo = ^H2/^i2co (expressed in cm ^ K ^ km ^ s), 
both near the GC (r 0) and either in the inner Galactic disk (r ~ 2 — 7 kpc) or in the vicinity of the Sun (r ^ tq). 



Aco(r~0) 


^co (r ~ 2 - 7 kpc) 


Xco{r ~ tq) 


Reference 


~ (2 - 7) X 10"'' 


~ 2.2 X 10^" 




Sodroski et al. (1995) 


~ 9 X 10^^ 




~ 3 X 10^" 


Arimoto et al. (1996) 


~ (2.4-7.2) X 10^^ 


~ 3 X 10^° 




Oka et al. (1998) 


~ (2 - 4) X 10^^ 




~ 1.5 X 10^° 


Strong et al. (2004) 



The use of the standard virial mass as a tracer of 
the actual mass of molecular clouds was called into ques- 
tion by Oka et al. (1998), who argued that molecular 
clouds near the GC are not gravitationally bound, but 
instead confined by the high external pressure from hot 
gas and magnetic fields. Using the correct expression 
for the virial mass of pressure-confined clouds, they re- 
calibrated the CO-luminosity-virial-mass relation for GC 
molecular clouds, whereupon they obtained Xco ^ 2.4 x 
10^^ cm-2 K"^ km-i s near the GC. This value of Xco 
is one order of magnitude lower than the value Xqq ~ 
3 X 10^° cm~^ km"-'^ s applying to gravitationally- 
bound molecular clouds in the Galactic disk. With their 
new value of Xco, Oka et al. (1998) estimated the H2 mass 
inside r ~ 375 pc at ~ 2 x 10^ Mq. However, they pointed 
out that their derived values were lower limits, which could 
underestimate the true conversion factor and H2 mass near 
the GC by up to a factor of 3. 

More recently. Strong et al. (2004) showed that an 
outward radial gradient in Xqo was required to rec- 
oncile the predictions of their cosmic-ray propagation 
code GALPROP with the 7-ray profiles measured by 
EGRET/COMPTON, when they adopted the molecular 
and atomic gas distributions inferred from CO and Hi 
surveys and assumed that the distribution of cosmic-ray 
sources follows the observed distribution of pulsars. They 
did not attempt to derive the precise r-dependence of Xqo 
implied by the measured 7-ray profiles, but they showed 
that a good fit to the EGRET data could be obtained 
with Xco dropping from ^ 1.5 x 10^° cm~^ km"-'^ s 
at r = r0 to ~ 4 X 10^^ cm^^ K^^ km~^ s at r = 2 kpc. 
Extrapolating their curve Xcoi''') in to the origin yields 
Xco ~ (2 - 4) X 10^^ cm-2 km^^ s at r (see 
their Figure 2). 

All the above estimates of the ^^CO-to-H2 conversion 
factor are summarized in Table [2l It is clear that the value 
of Xco near the GC remains fairly uncertain. More im- 
portantly, the very notion that Nn^ is proportional to 
W^i2co may be questionable (e.g., Dahmen et al. 1998). 
Nevertheless, the fact that three completely different meth- 
ods to determine Xco near the GC (far-infrared emission 
from dust, virial masses of molecular clouds and 7-ray emis- 
sion from cosmic rays) broadly converge to a common range 
of values lends some credence to the presumed proportion- 
ality between and W^i2co- For the following, we adopt 
Xco('' ~ 0) = 5 X 10^^ cm~^ km~^ s as a reference 
value, which we consider uncertain by a factor ~ 2. 

Once the value of Xco is known, the H2 spatial distribu- 
tion can in principle be inferred from CO emission measure- 
ments. The overall disposition of CO emission in the Galaxy 
was brought to light by the early large-scale Galactic CO 
surveys, which revealed a sharp peak inside a few 100 pc 
from the GC, a deep depression between ^ 1.5 kpc and 
3 kpc and a ring of enhanced emission (now designated the 



Galactic disk molecular ring) between ~ 3.5 kpc and 6 kpc 
(see Combes 1991 for a review). 

Sanders et al. (1984) presented a full-coverage Galactic 
^^CO (J = 1 — > 0) survey, on the basis of which they 
constructed an axisymmetric model of the H2 distribu- 
tion in the Galaxy. According to their model (rescaled to 
tq = 8.5 kpc), molecular gas in the GB is confined to a 
disk with radius ~ 1.3 kpc and nearly constant thickness 
~ 70 pc, seen edge-on from the Sun and tilted counter- 
clockwise by ~ 7° with respect to the Galactic planeQ (so 
that its midplane passes from negative latitudes in the first 
quadrant to positive latitudes in the fourth quadrant). This 
tilted disk is prolonged by a rarefied molecular layer with 
approximately the same thickness, which smoothly joins up 
with the Galactic disk molecular ring. Taking a constant 
conversion factor Xco = 3.6 x 10^° cm~^ km~^ s for 
the whole Galaxy, Sanders et al. (1984) found that the H2 
space-averaged density could be approximated by 



(nH2)(r) = (nH2)o(r) exp 



H 



(4) 



where (after rescaling to tq — 8.5 kpc) H ~ 
(70 pc)/(2 V2 In 2) ~ 30 pc out to r ~ 3 kpc, zo(r,6') 
is the local Galactic height of the H2 midplane (given by 
the tilted disk geometry inside 1.3 kpc) and (7^-112) o(^) is 
the H2 space-averaged density at zq, such that (?^H2)o(0 ~* 
450 pc) ~ 74 cm-^, (nH2)o(450 pc 1.3 kpc) ~ 12 cm'^ 
and (nH2)o(l-3 3 kpc) ~ 1.1 cm~^. The corresponding 
H2 masses inside r = 450 pc, between 450 pc and 1.3 kpc, 
and in the depression zone between 1.3 kpc and 3 kpc are 
listed in Tabled 

Burton & Liszt (1992) presented higher-resolution 
^■^CO ( J = 1 — > 0) observations toward the GC and pointed 
out that the central parts of the gaseous GB disk are in 
fact closely aligned with the Galactic plane. More specif- 
ically, they showed that the observed {l,b,v) pattern of 
^^CO emission in the GB could be understood in the frame- 
work of a tilted and warped H2 disk model with the fol- 
lowing characteristics (rescaled to = 8.5 kpc): Out to 
r ~ 170 pc, the H2 disk is flat, parallel to the Galactic 
plane and ~ 30 pc thick (Gaussian scale height ~ 13 pc). 
Between ~ 170 pc and 1.5 kpc, the H2 disk flares linearly 
to a thickness ~ 300 pc (uncertain value), and this flar- 
ing is accompanied by a warp such that the H2 midplane 
becomes tilted with respect to the Galactic plane, by an an- 
gle that varies sinusoidally with Galactic azimuth between 
— 13° at 9 = 45° (so that < 0° in the near, flrst-quadrant 
sector) and -^-13° at 6* = 225° (6 > 0° in the far, fourth- 
quadrant sector) (see their Figure 5). Overall, the warped 



Here, and in the rest of the paper, the Galactic plane is 
defined as the plane 6 = 0°, which is slightly different from the 
plane of the solar circle. 
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Region of space 



Radial range 



Ho mass 



Reference 



Tilted GB disk : inner part 
outer part 

Depression zone 
GC molecular disk 
"Expanding molecular ring" 
GB region -1?5 <l<3°5 
GB region K| < 2?5 
Nuclear bulge 



- 450 pc 

450 pc - 1.3 kpc 

1.3 — 3 kpc 

~ - 150 pc 

straddling r ~ 180 pc 

~ — 375 pc (off-centered) 

~ - 375 pc 

~ - 280 pc (off-centered) 



~ 1.75 X lO"* Mq 
~ 2.1 X 10* Mq 
~ 9 X 10^ Mq 
~ 2.4 X lO'' Mq 
~ 5 X 10** Mq 
~ (2 - 6) X lO'' Mq 
~ (2 - 6) X lO'^ Mq 
Dense gas : ~ (1.2 - 
Thin gas : ~ (0.7 - 



- 6.4) X 10^ Mq 
1.4) X 10^ Mq 



Sanders et al. (1984) 



Sofue (1995a) 
Sofue (1995b) 
Sodroski et al. (1995) 
Oka et al. (1998) 
Dahmen et al. (1998) 



H2 disk appears tilted out of the Galactic plane and in- 
clined to the line of sight, very much like the tilted Hi disk 
of Burton & Liszt (1978) (see section 1^75)) . The warped H2 
disk also bears a resemblance to Sanders et al.'s (1984) 
tilted H2 disk, but it differs by its substantial flaring and 
by the direction of maximum tilt {9 — 45° - 225°, as op- 
posed to 9 — 90° - 270° in Sanders et al.). Interestingly, 
the H2 space-averaged density in the warped disk can be 
expressed by Equation^ with ~ 13 pc inside r ~ 170 pc 
and H ~ 0.088 {r - 22 pc) between ~ 170 pc and 1.5 kpc, 
and with (nH2)o oc 1/H (such that the H2 column density 
through the disk has a constant, albeit unspecified, value). 
Besides its geometric parameters, the warped H2 disk model 
also possesses two kinematic parameters describing rotation 
and expansion motions, respectively. 

Let us now focus on the innermost Galactic regions, in- 
side a few 100 pc from the GC, where molecular gas tends to 
concentrate. Sofue (1995a) investigated the H2 morphology 
and kinematics inside |^| ~ 1° {r± ~ 150 pc), based on the 
AT&T BeU Laboratories survey of ^^CO ( J = 1 ^ 0) emis- 
sion conducted by Bally et al. (1987). He adopted Xqo = 
9.2 X 10^^ cm~^ K"-'^ km~^ s (from a slightly modified ver- 
sion of Equation [S]) together with Wi2co/M^i3co = 6.2 
(measured intensity ratio averaged over \l\ < 20°). The 
most salient continuous features in the longitude-velocity 
{I, v) diagrams of the GC molecular disk are two ap- 
parently rigidly-rotating ridges. Sofue (1995a) interpreted 
these ridges as the {I, v) traces of two dense material arms, 
which together would form a rotating ring with mean ra- 
dius ~ 120 pc and vertical thickness ~ 15 pc (except for 
the massive molecular complexes around Sgr B and C, 
whose thickness reaches ^ 30 — 60 pc) . An analysis of the 
velocity-integrated intensity maps yields a total H2 mass 
- 2.4 X 10'^ Mq for the disk, with - 1.1 x 10^ Mq in 
Arm I and ~ 0.8 x 10^ Mq in Arm II, and an H2 mass 
~ 2.8 X 10^ Mq for the entire region r± < 150 pc. By ex- 
amining a somewhat wider section of Bally et al.'s (1987) 
survey, Sofue (1995b) estimated that the "expanding molec- 
ular ring" that surrounds the GC molecular disk has a mean 
radius ~ 180 pc (as already established before), a radial 
thickness ^ 15 pc, a vertical thickness ~ 100 pc and an H2 
mass ^ 5 X 10^ Mq (80% of which belongs to the region 
rj_ < 150 pc). 

Dahmen et al. (1998) gained additional information on 
the spatial structure and physical state of molecular gas 
near the GC by analyzing a SMWT C^^O (J = 1 ^ 0) 
survey of the central region — 1?05 < I < 3?6, in conjunc- 
tion with the SMWT ^^CO (J = 1 ^ 0) survey of the 
broader region —12° < / < 13° presented by Bitran et 



al. (1997). For both isotopomers|3 they performed radia- 
tive transfer calculations in the Large Velocity Gradient 
(LVG) approximation, and they adopted abundance ratios 
12CO/H2 = 10-4 and CISO/H2 = 4 X 10"^. Their calcu- 
lations enabled them to estimate the kinetic temperature, 
volume density and total mass of molecular gas in the nu- 
clear bulge, which they assumed to extend over the longi- 
tude range -1?5 < / < 2?25 (projected half-size ~ 280 pc). 
They found that most of the C^^O emission is likely to orig- 
inate in molecular gas with kinetic temperature T ~ 50 K 
and H2 density ~ lO"^-^ cm~'^, and they argued that 
the associated H2 mass in the nuclear bulge is comprised 
between ~ 1.2 x 10^ Mq and 6.4 x lO'^ Mq. Furthermore, 
they noted significant differences between the C^^O and 
"'^^CO intensity maps, which they interpreted as evidence 
for the presence of a widespread, high-temperature and low- 
density gas component, detected in ^^CO but not in C^^O. 
This thin gas would have T ~ 150 K, riHa ~ 10^'^ cm~'^ and 
an H2 mass in the nuclear bulge between ~ 0.7 x 10'' Mq 
and twice that value. 

Another powerful method to probe the physical condi- 
tions in molecular gas near the GC is to compare emission 
lines from different J levels of ^^CO. In this spirit. Oka et 
al. (1998) mapped the region \l\ < 2?5 in ^^CO ( J = 2 ^ 1) 
emission with the NRO survey telescope, and they resorted 
to previous ^^CO (J = 1 — > 0) maps from Bitran (1987) 
for comparison. Their analysis of the ^^CO (J = 2 1) 
/ ^^CO ( J = 1 ^ 0) intensity ratio indicates that the CO 
luminosity of the innermost ~ 400 pc is dominated by emis- 
sion from low-density gas with riHa ~ 10^'^ cm~'^. To recon- 
cile this finding with the inference from CS, NH3 and HCN 
observations that high-density gas with riH2 ^ cm"'^ 
prevails. Oka et al. (1998) advanced the view that molec- 
ular gas actually exists in two distinct components: a low- 
density (tt-Hj ^ 10'^ cm~^) component present in "diffuse" 
clouds with a large filling factor, and a high-density (nua ^ 
10^ cm~'^) component present in "clumps" with a small 
filling factor. This view of a two-component molecular gas 
in the GC region has found support in several subsequent 
studies relying on other tracer molecules (e.g., Rodriguez- 
Fernandez et al. 2001; Oka et al. 2005; Magnani et al. 2006). 
Let us simply mention that the warm (T 150 K), and 
presumably low-density, component detected by Rodriguez- 
Fernandez et al. (2001) would on average represent ~ 30% 
of the total molecular gas. 



Isomers having the same number of isotopically distinct 
atoms but differing in the positions of these atoms. The term 
"isotopomer" results from the contraction of the words "isotope" 
and "isomer". 



6 



K. Ferriere et al.: Spatial distribution of interstellar gas in the Galactic bulge 



Table 4. Estimated values of the mass of interstellar H2 in specific regions of the GB, deduced from ^^CO measurements 
(or from ^'^CO measurements combined with an assumed ^^CO-to-^'^CO intensity ratio) together with Xqo = 5 x 
IQi^ cm-2 K~i km~i s. 



Region of space 


Radial range 


H2 mass 




Original reference" 


Tilted GB disk : inner part 


- 450 pc 


~ 2.4 X 10' 


Mq 


Sanders et al. (1984) 


outer part 


450 pc - 1.3 kpc 


~ 2.9 X 10^ 


Me 




Depression zone 


1.3 — 3 kpc 


~ 1.3 X lO'' 


Mq 


): 


GC molecular disk 


~ - 150 pc 


~ 1.3 X lO'^ 


Mq 


Sofue (1995a) 


"Expanding molecular ring" 


straddling r ~ 180 pc 


~ 2.7 X 10^ 


Mq 


Sofue (1995b) 


GB region -1?5 < Z < 3?5 


~ — 375 pc (off-centered) 


~ 4.4 X 10'' 


Mq 


Sodroski et al. (1995) 


GB region \l\<2°5 


~ - 375 pc 


~ 4.2 X 10^ 


Mq 


Oka et al. (1998) 


" Reference for the eorrespondinj 


; mass estimate given in Table [S] which we rescaled to our adopted value of Xqq. 



Higher-J transitions of ^^CO also provide a valuable 
diagnostic tool for molecular gas. For instance, Martin et 
al. (2004) surveyed the region -1?3 < I < 2° in ^^CO (J = 
4^3) and ^^CO (J = 7 ^ 6) emission with the AST/RO 
telescope, and they analyzed their ^^CO ( J = 4 — + 3) and 
12C0 (J = 7 6) maps together with existing ^^CO (J = 
1 0) and i^CO (J = 1 ^ 0) maps from Stark et al. 
(1988) and Bally et al. (1987, 1988). Like Dahmen et al. 
(1998), they worked with the LVG approximation and they 
took ^2CO/H2 = ID"**. What emerges from their study is 
that the molecular gas kinetic temperature decreases from 
T > 70 K at the edges of cloud complexes to T < 50 K in 
their interiors, and that the H2 density spans the validity 



range of their LVG analysis (~ 10 



2.5 



10 



,4.5 



cm 



typical values 71^2 ~ 10 cm inside cloud complexes. 

Finally, Sawada et al. (2004) devised a purely obser- 
vational method to deconvolve sky maps of CO emission 
into a face-on view of the molecular gas distribution near 
the GC, based on a quantitative comparison between the 
2.6-mm ^^CO emission line and the 18-cm OH absorption 
line. Their face-on map exhibits a strong central condensa- 
tion (corresponding to the CMZ), which is elongated along 
an axis inclined by ^ 70° to the line of sight (so that its 
near end hes at positive longitudes), ~ 500 pc x 200 pc in 
size, and possibly composed of a pair of arms (as suggested 
by Sofue 1995a). The so-called "expanding molecular ring" 
lies at the periphery of the central condensation, without 
clearly separating from it. It, too, is elongated and inclined 
(by an uncertain angle < 70°) toward positive longitudes, 
though it does not actually stand out as a coherent entity 
in the face-on map. Instead, Sawada et al.'s (2004) results 
tend to substantiate the idea that the "expanding molecular 
ring" arises from the projection of interstellar gas moving 
along highly elongated orbits in response to the gravita- 
tional potential of the Galactic bar (see section[3]for further 
details) . 

Table [3] provides a summary of all the H2 masses esti- 
mated in this subsection. Clearly, the H2 mass derived by 
Sanders et al. (1984) for the innermost region is significantly 
greater than the other estimates. This is only because these 
authors used a much larger value of the conversion factor, 
Xqo, more typical of the Galactic disk at large. If we now 
rescale all the H2 masses deduced from ^^CO measurements 
(or from ^^CO measurements combined with an assumed 
i^CO-to-i^CO intensity ratio) in Table [3] to our adopted 
reference value Xqo = 5 x 10^^ cm^^ K^^ km^^ s, we ob- 
tain the H2 masses listed in Table ID These masses agree 
with each other to within a factor of ^ 2. Moreover, they 
grossly match the total hydrogen masses inferred from in- 



frared dust measurements by Launhardt et al. (2002) (see 
Tabled]), consistent with the notion that most of the inter- 
stellar gas in the CMZ is in molecular form. 

2.3. Atomic gas 

Virtually all our observational knowledge on the interstel- 
lar atomic gas near the GC stems from spectral studies of 
the Hi 21-cm line in emission or in absorption. The early 
NRAO surveys of Hi 21-cm emission over the longitude 
range -11° < / < 13° led Burton & Liszt (1978) and Liszt 
& Burton (1980) to propose a tilted disk model for the Hi 
spatial distribution inside r ~ 2 kpc. The model parameters 
), with were determined iteratively by generating synthetic 21-cm 
emission spectra and adjusting them to the observed spec- 
tra. 

In Burton & Liszt's (1978) model (rescaled to tq = 
8.5 kpc), the Hi disk is axisymmetric, 1.3 kpc in radius, 
200 pc thick (Gaussian scale height of 85 pc), tilted by 
22° out of the Galactic plane (with the first/fourth quad- 
rant side at negative/positive latitudes) and inclined by 78° 
to the plane of the sky (with the near/far side at nega- 
tive/positive latitudes), such that the normal to the Hi disk 
forms a 25° angle with the z-axis. The Hi space-averaged 
density in the disk is given by a Gaussian function of dis- 
tance to the Hi midplane, Zd, 



("-Hi)(zd) = (0.39 cm ^) exp 



Zd 



85 pc 



(5) 



which implies a total Hi mass in the disk ~ 1.1 x 10^ Mq. 
Finally, the Hi gas within the disk has rotation and expan- 
sion motions of comparable magnitudes. 

Liszt & Burton (1980) argued that the observed kine- 
matics of the Hi gas could be more plausibly explained by 
a tilted elliptical disk model, in which the gas moves along 
closed elliptical orbits. Thus, in Liszt & Burton's (1980) 
model (again rescaled to =8.5 kpc), the Hi disk is ellip- 
tical, with semi-major axis 1.6 kpc and axis ratio 3.1 : 1. Its 
thickness remains 200 pc, but the tilt and inclination angles 
are now 13?5 to the Galactic plane and 70° to the plane of 
the sky (entailing a 24° angle between the disk normal and 
the z-axis), and the major axis is oriented at 48? 5 clockwise 
to the "inclined a;- axis" (corresponding to a 51? 5 angle to 
the hne of sight) (see their Figure 1). The Hi space-averaged 
density in the disk is still given by Equation[5l but the total 
Hi mass in the disk now amounts to ~ 5.2 x 10^ AIq. 

A different modification to Burton & Liszt's (1978) orig- 
inal model can be found in the paper by Burton & Liszt 
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Region of space 



Radial range 



Hi mass 



Reference 



Tilted GB disk : 
Axisymmetric 
Elliptical 
Axisymmetric 
Inner part of GB disk 
180-pc molecular ring 



- 1.3 kpc 

- 1.6 kpc / - 0.52 kpc 

See Equation |6] 

~ — 300 pc (off-centered) 

straddling r ~ 175 pc (off-centered) 



' 1.1 X lO'^ Mq 

■ 5.2 X 10'' Mq 

■ 2.2 X 10® Mq 

■ 2.3 X lO'' Mq 
' 8 X 10^ Mq 



Burton & Liszt (1978) 
Liszt & Burton (1980) 
Burton & Liszt (1993) 
Rohlfs & Braunsfurth (1982) 



(1993). There, the Hi spatial distribution within the GB 
disk is modeled in a slightly more complicated manner, with 
a Hi space-averaged density depending on both r and Zd'. 



{nm){r., Zd) 



(0.82 cm"^) exp 



-h(0.10 cm"^) exp 



X exp 



1 



Zd 



2 V50 pc 




(6) 



The corresponding FWHM thickness is ~ 118 pc and the 
total Hi mass is ~ 2.2 x 10^ Mq. However, the emphasis 
of the paper is placed on the kinematic aspects of atomic 
gas in the GB, and as the authors themselves recognize. 
Equation [6] may have little intrinsic significance. Therefore, 
we do not regard the associated mass estimate as particu- 
larly trustworthy. 

Rohlfs & Braunsfurth (1982) presented a more refined 
description of the morphology and kinematics of atomic gas 
in the central parts of the Hi GB disk. They analyzed 21- 
cm emission as well as absorption spectra obtained with the 
Effelsberg telescope at longitudes |/| < 1?5, and they made 
comparisons between both types of spectra in order to lo- 
cate the major Hi features relative to the GC. Like Burton 
& Liszt (1978), they started from the premise that the ve- 
locity field includes rotation and expansion. Furthermore, 
they used an axisymmetric rotation curve inferred from 
a Galactic mass model as an independent distance indi- 
cator, and they assumed that the major Hi features are 
nearly circular (though not necessarily concentric with the 
GC). Because of these stringent assumptions, largely re- 
futed by more recent work, too much reality should not 
be ascribed to the precise quantitative estimates provided 
by Rohlfs & Braunsfurth (1982). Nevertheless, a couple of 
important points deserve mentioning. First, although the 
Hi GB disk as a whole is markedly tilted, its central parts 
are approximately parallel to the Galactic plane (as is the 
case for the warped H2 disk of Burton & Liszt 1992; see 
section [521) . Second, the 180-pc molecular ring is visible 
in 21-cm emission, without however being prominent. For 
reference, Rohlfs & Braunsfurth (1982) estimated that the 
inner part of the Hi GB disk has a radius ^ 300 pc, a 
thickness ~ 90 pc and an Hi mass ^ 2.3 x 10^ Mq, while 
the 180-pc molecular ring has a mean radius ~ 175 pc, a 
vertical thickness ~ 60 pc and an Hi mass ~ 8 x 10^ Mq. 

Unfortunately, the atomic gas does not benefit from the 
wide variety of diagnostic spectral lines or from the vast 
number of high-resolution surveys as are presently avail- 
able for the molecular gas. Further genuine progress in 
our understanding of interstellar atomic gas near the GC 



will probably have to await the release of new data from 
the GC extension of the Southern Galactic Plane Survey 
(McClure-Griffiths, private communication). In the mean- 
time, one should bear in mind that the spatial distribution 
and kinematics of the atomic gas might not be as drastically 
different from those of the molecular gas as traditionally be- 
lieved. Liszt & Burton (1996) advanced several cogent argu- 
ments pointing to the congruence of both gases in position 
and in velocity, and they attributed their apparent dispar- 
ities to differences in the observational conditions (angular 
resolution, sky coverage and signal-to- noise ratio) of Hi ver- 
sus CO, in the thin versus thick optical-depth regime and in 
the presence versus absence of absorption (see also Liszt & 
Burton 1978). In addition, they emphasized that the kine- 
matic features appearing in the Hi and CO spectra (notably 
the "expanding molecular ring") do not correspond to dis- 
crete material bodies, but instead result from kinematic 
projection effects. 

Again, we close up the subsection with a summary table 
(Table [5]) including all our mass estimates. For complete- 
ness, we can add that the Hi mass inside r ~ 300 pc is 
~ 5.8 X 10^ Mq in Burton & Liszt's (1978) model and 
~ 6.6 X 10^ Mq in their (1993) model, as opposed to 
~ 3.1 X 10^ Mq in Rohlfs & Braunsfurth (1982). Hence, 
there is roughly a factor of 5 discrepancy between Burton 
& Liszt on the one hand and Rohlfs & Braunsfurth on the 
other hand. The same factor exists, for the entire Hi GB 
disk, between the old and new models of Burton & Liszt, 
but as we already mentioned earlier, the Hi mass deduced 
from their 1993 model should not be taken too seriously. If 
we disregard this mass estimate, we may conclude from a 
comparison between Tables |4] and [5] that the Hi mass of the 
entire GB disk represents but ~ 10% — 20% of its H2 mass. 



2.4. Ionized gas 

Radio signals from pulsars and other (Galactic and extra- 
galactic) compact sources provide a unique source of infor- 
mation on the ionized component of the interstellar gas. 
Cordcs & Lazio (2002) assembled all the useful data (dis- 
persion measures, scattering measures and independent dis- 
tance estimates) available at the end of 2001 to construct 
a new non-axisymmetric model of the spatial distribution 
of interstellar free electrons in the Galaxy. This "NE2001 
model" , which incorporates separate multi-wavelength data 
on the spiral structure of the Galaxy, the local interstellar 
medium and the GC region, supersedes the earlier mod- 
els of Cordes et al. (1991) and Taylor & Cordes (1993). 
It includes a smooth large-scale component, which consists 
of two axisymmetric disks and five spiral arms, a smooth 
GC component, a contribution from the local interstellar 
medium, and individual clumps and voids. 
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Retaining solely the smooth components present inside 
3 kpc, we can write the free-electron space-averaged density 
in the GB as 

(no)(r) (nc)i(r, z) + (nc)2(r, z) + (nc)3(r-, z) , (7) 
where 



(ne)i(r,z) = (0.05 cm-3) 



cos vr 



2Li 



u{Li - r) 



X sech 



(8) 



with Li = 17 kpc and Hi = 950 pc, for the outer thick disk 
(u is the unit step function), 



)2(r,z) = (0.09 cm-3) exp 



r — L2 
L2/2 



X sech I — 
H2 



(9) 



with L2 — 3.7 kpc and H2 = 140 pc, for the inner thin 
disljfl, and 



(nc)3(r, z) = (10 cm ^) exp 



x'^ + iy- Vz) 



21 



7^2 
^3 



X exp 



- ^3) 



21 



^1 



(10) 



with 2/3 = —10 pc, Z3 = — 20pc, L3 = 145 pc and Hj, = 
26 pc, for the GC component [j Strictly speaking, the offset 
ellipsoid exponential in Equation [TO] should be truncated 
to zero for arguments smaller than —1. However, such a 
discontinuity was introduced into the model to reflect the 
abrupt changes observed in the distribution of scattering 
diameters for OH masers in the GC rather than an abrupt 
drop-off in the free-electron density. Since we are interested 
only in the free-electron density (as opposed to both the 
density and its fluctuations), we are entitled to relax the 
constraint arising from scattering diameters and omit to 
truncate the ellipsoid exponential in Equation [TO] (Lazio, 
private communication). Then the GC component remains 
dominant near the midplane out to r ~ 335 pc. The rest of 
the GB is largely dominated by the thick disk, except for a 
narrow wedge starting at (r, z) ~ (2.3 kpc, 0) and opening 
outward, wherein the thin disk takes over. 

To convert the above free-electron density into an 
ionized-hydrogen density, we follow conventional wisdom 
and assume that the molecular and atomic media (discussed 
in sections 12.21 and 12.31 respectively) can be considered as 
completely neutral, while the ionized medium can be di- 
vided into a warm ionized medium (WIM), where hydro- 
gen is entirely ionized and helium completely neutral, and 
a hot ionized medium (HIM), where hydrogen and helium 
are both fully ionized. If we further denote by /him the 
fraction of ionized gas belonging to the HIM and remember 
that helium represents 10% by number of hydrogen (see 
beginning of section [2]), we find that the ionized-hydrogen 



^ The expression of g2{r) given in Table 2 of Cordes & Lazio 
(2002) contains a typo: Aa should be A2 in the numerator and 
A2/2 in the denominator (Lazio, private communication). 

® The {x, y, z) coordinates used by Cordes & Lazio (2002) cor- 
respond to (y,x,z) in our coordinate system (see Figure[l|. 



space-averaged density is related to the free-electron space- 
averaged density through 



(n-H+ 



1 



1 + 0.2 /] 



HIM 



(11) 



The HIM fraction, /him, is quite uncertain, but fortunately, 
its exact value has little impact on the inferred H+ density, 
(nH+), which varies between (rio) if /him — and 0.83 (ric) 
if /him = 1- 

The total H+ mass in the interstellar GB and the con- 
tributions from the thick disk, thin disk and GC compo- 
nent are obtained by integrating Equations [5] [5] and [TOl 
divided by the factor (1 + 0.2 /him), out to r = 3 kpc. For 
reference, if /him = 0, the thick disk has an H+ mass of 
6.4 X 10^ Mq inside 3 kpc, with 1.8 x 10^ Mq inside 500 pc 
and 1.6 x 10^ Mq inside 150 pc; the thin disk has an H+ 
mass of 8.3 X 10^ Mq inside 3 kpc, with 1.8 x 10"* Mq inside 
500 pc and 1.0 x 10^ Mq inside 150 pc; and the GC compo- 
nent has a total H+ mass of 7.5 x 10^ Mq , with 5.0 x 10^ Mg 
inside 150 pc and negligible amounts outside 500 pc. The 
resulting total H+ masses inside 150 pc, 500 pc and 3 kpc 
are displayed in Table [S] A comparison with our previous 
tables (Table [T] or Table |4] with an assumed Hl-to-H2 mass 
ratio ~ 10% — 20%) suggests that ionized gas accounts for 
only ~ 4% — 5% of the total interstellar gas in the region 
r < 150 pc, ~ 4% — 8% in the region r < 500 pc, and 
^ 30% — 50% on average over the entire GB (r < 3 kpc). 
Not surprisingly, the fraction of ionized gas increases not 
only with |z| (as indicated by the comparatively large scale 
height in Equation [8]) , but also with r. 

Lazio and Cordes (1998) discussed the nature of the 
medium responsible for the observed scattering. They pro- 
posed that scattering arises in thin layers on the surfaces 
of molecular clouds. These layers would either be the pho- 
toionized skins of molecular clouds, with Te ~ lO"' K and 
TT-e ^ 10'^ cm~^, or the interfaces between molecular clouds 
and the hot (- 10^ K) ambient gas, with ~ (lO^^-lO^) K 
and Tie ~ (5 — 50) cm^'^. In both cases, the scattering 
medium would have a small filling factor. However, the 
scattering medium possibly constitutes but a small part of 
the ionized medium; therefore, its morphology and physical 
parameters are not necessarily representative of the ionized 
medium in general. 

The best source of information on the morphology and 
physical parameters of the warm component of the ionized 
medium near the GC comes from the thermal (free-free) 
component of the radio continuum emission and from radio 
recombination lines. In general, thermal continuum mea- 
surements are used to image the ionized gas and to estimate 
its electron density and temperature, while recombination 
line measurements make it possible to investigate the kine- 
matics of the ionized gas and also to estimate its electron 
temperature. 

As explained in the review by Mezger & Pauls (1979), 
the thermal radio continuum emission is produced in an 
extended ionized medium (here referred to as the WIM) 
and in individual Hii regions. On the sky, the thermal ra- 
dio emission appears closely correlated with far-infrared 
dust emission, and less well correlated with CO emission, 
which is asymmetric (see section 12.21) , more extended in 
longitude and less extended in latitude. To quote Mezger 
& Pauls (1979), "this suggests that the [WIM] may be 
ionization bounded along the Galactic plane, but density 
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Table 6. Estimated values of the mass of interstellar H+ (total mass or partial contribution from the WIM or 
HIM/VHIM) in specific regions of the GB. The total H+ masses deduced from the NE2001 model for the free-electron 
density (Cordes & Lazio 2002) are those obtained under the extreme assumption that 100% of the ionized medium is 
warm (so that (71h+) ~ ("-c))- If, in reality, a fraction /him of the ionized medium is hot, these masses should be divided 
by the factor (1 + 0.2 /him) (see Equation [TT|) . 



Region of space 


Radial range 


Total 


mass 
WIM 


HIM/VHIM 


Reference 


Galactic bulge : innermost 


- 150 pc 


~ 6.6 X lO'' Mq 






Cordes & Lazio (2002) 


intermediate 


- 500 pc 


~ 2.6 X 10" Mq 








total 


- 3 kpc 


~ 7.3 X lO'^ Mq 








Smaller ellipsoid 


~ - 47 pc 




~ 1.2 X lO'^ Mq 




Mezger & Pauls (1979) 


Larger ellipsoid 


~ 0-112 pc 




~ 4.7 X lO'^ Mq 






Galactic bulge 


- 3 kpc 






~ 8 X 10® Mq 
~ 1.2 X lO'^ Mq 


Snowden et al. (1997) 
Almy et al. (2000) 


10-keV plasma : innermost 


- 150 pc 






~ 5.5 X 10'' Mq 


j- Yamauchi et al. (1990) 
1 Koyama et al. (1996) 


total 


~ - 250 pc 






~ 9.5 X 10'' Mq 


55 



bounded perpendicular to it." As a first rough approxi- 
mation, Mezger & Pauls (1979) modeled the WIM by the 
superposition of two oblate ellipsoids with the following 
characteristicfQ (rescaled to tq — 8.5 kpc): the larger el- 
lipsoid is ~ (225 pc)^ X 90 pc in size and has Te — 5000 K, 
rig ~ 8 cm~^ and Mh+ ~ 4.7 x 10^ Mq, while the smaller 
ellipsoid is ~ (95 pc)^ x 55 pc in size and has Te ~ 5000 K, 
rie ~ 18 cm-3 and Mh+ ~ 1.2 x lO'^ Mq (see Table E]). 
Altogether, the central electron density is ~ 26 cm~^ and 
the total H+ mass is ~ 5.9 x 10^ Mq. For comparison, in the 
NE2001 model of Cordes & Lazio (2002), the GC compo- 
nent has an ellipsoid exponential distribution with FWHM 
size ~ (241 pc)^ x 43 pc (comparable in radial extent, but 
significantly more oblate than the larger ellipsoid of Mezger 
& Pauls 1979), a central electron density of 10 cm~'^ and a 
total H+ mass of (7.5 x 10^ Mq)/{1 + 0.2 /him)- Evidently, 
the good agreement found for the H''" masses is a little for- 
tuitous. 

A number of studies on the WIM near the GC have been 
carried out since the review of Mezger & Pauls (1979). For 
instance, Mehringer et al. (1992, 1993) observed the Sgr B 
region with the VLA, both at several wavelengths in the ra- 
dio continuum and in the HllOa radio recombination line. 
They estimated the r.m.s. electron density outside compact 
Hii regions at ~ 80 cm~^ in Sgr Bl and ~ 60 cm~'^ in Sgr 
B2. 

Besides the WIM, the ionized medium near the GC con- 
tains a hot component (referred to as the HIM) which can 
be detected through its X-ray thermal emission. Maps of 
the diffuse X-ray background in the (0.5 — 2.0) keV energy 
band from the ROSAT all-sky survey reveal an extended 
zone of enhanced emission in the general direction of the 
GC (Snowden et al. 1997). Although part of the enhance- 
ment can be attributed to the nearby Loop I SB, Snowden 
et al. (1997) argued, based on the latitude profile of the ex- 
cess emission, on the absorption trough running along the 



In the rest of this section, all density estimates deduced 
from emission measures implicitly assume that the phase under 
consideration occupies all the interstellar volume. If the consid- 
ered phase actually has a volume filling factor 0, each of our 
density estimates corresponds in fact to the geometric mean 
between the true and space-averaged densities, and the actual 
true/space-averaged density is obtained by dividing/multiplying 
our density estimate by y^. 



Galactic plane and on the deep shadows cast by relatively 
distant molecular clouds, that the bulk of the enhancement 
arises from a bulge of hot, X-ray emitting gas around the 
GC. They went on to propose a crude model for the X-ray 
bulge, which relies on the (0.5 - 2.0) keV ROSAT data cor- 
rected for foreground absorption by interstellar Hi and H2, 
and which assumes that the X-ray emitting gas is in col- 
lisional ionization equilibrium. In their model, the hot gas 
has a temperature of 4 x 10^ K, an electron space-averaged 
density (again rescaled to Tq = 8.5 kpc) 



("c)him(^-2) = (0.0034 cm-3) u(6kpc-r) 



X exp 



2 kpc 



(12) 



and hence an H+ mass ~ 8 x 10® Mq inside r — 3 kpc. 

As follow-up work, Almy et al. (2000) analyzed the 
shadows cast by a distant molecular cloud complex against 
the X-ray enhancement in the | keV and 1.5 keV ROSAT 
bands, and came to the conclusion that a significant frac- 
tion of the observed radiation in that direction is indeed 
emitted beyond the cloud complex, most likely in an X-ray 
bulge surrounding the GC. They also derived a simple, al- 
beit more realistic, model for the X-ray bulge, in which the 
X-ray emitting gas is an adiabatic polytrope (P oc p^^"^) in 
hydrostatic balance (in the Galactic gravitational potential 
of Wolfire et al. 1995). Like Snowden et al. (1997), they as- 
sumed coUisional ionization equilibrium and they accounted 
for foreground absorption. They found that the best fit to 
the I keV ROSAT data is reached for a central gas tem- 
perature of 8.2 X 10® K and a central electron density of 
0.011 cm-3. 

Almy et al. (2000) provided no analytical expressions 
to describe the spatial variation of the hot gas parameters. 
Nonetheless, all the useful expressions can be retrieved by 
solving the hydrostatic equation, — VP— p V(/> = 0, together 
with the polytropic assumption, P = K p^l^^ and Wolfire 
et al.'s (1995) representation of the Galactic gravitational 



10 



K. Ferriere et al.: Spatial distribution of interstellar gas in the Galactic bulge 



potential, 
0(r,z) = 



-(225 km s"^)^ 




where Ci = 8.887 kpc, ai = 6.5 kpc, &i = 0.26 kpc, C2 = 
3.0 kpc, 02 = 0.70 kpc, C3 = 0.325, 03 = 12 kpc and 
r/i = 210 kpc. In this manner, it is found that the electron 
space-averaged density of the HIM can be written in the 
form 

K)him(^,z) = |(0.011 cni-'^f/'^ - (1.73 X 10"" cm"'^ s^) 
X [</.(r,z)- 0(0,0)]}''' (14) 



and that the 11+ mass inside r = 3 kpc (with |z| < 5 kpc) 
amounts to ~ 1.2 x 10'' Mq. 

Hence, the mass of hot H+ predicted by Almy et al.'s 
(2000) model for the region r < 3 kpc is 1.5 times that 
predicted by Snowden et al.'s (1997) cruder model. In view 
of Table [H] and the discussion below Equation [11] Almy 
et al.'s model suggests that hot gas globally represents ~ 
17% of all the ionized gas, and thus ^ 5% — 8.5% of the 
total interstellar gas, in the GB. However, little can be said 
about the spatial variation of the hot gas fraction, because 
the assumption of an adiabatic polytrope in hydrostatic 
balance is made for modeling convenience rather than on 
observational grounds. 

The ionized medium near the GC also contains an even 
higher-temperature component (sometimes called the very 
hot ionized medium or VHIM), revealed by its hard X- 
ray thermal emission, both in a free-free continuum and in 
characteristic lines from highly ionized elements. Koyama 
et al. (1989) reported the detection with the Ginga satel- 
lite of intense 6.7-keV Ka line emission from helium- like 
iron, the most abundant ionization state at a temperature 
of several keV. From the shape of the (supposedly free- free) 
continuum in the (2 — 18) keV energy band, they indepen- 
dently estimated the plasma temperature at kT ^ 10 keV 
(i.e., T ~ 10^ K). Yamauchi et al. (1990) then performed 
a two-dimensional mapping of the very hot 6.7-keV emit- 
ting plasma (outside the Galactic ridge). They found that 
its two-dimensional distribution on the plane of the sky 
could be fitted by an elliptical Gaussian with FWHM size 
~ 270 pc X 150 pc and major axis tilted clockwise by ~ 21° 
with respect to the Galactic plane. 

Subsequent, higher-energy- resolution observations in 
the (2—10) keV band with the ASCA satellite confirmed the 
presence of very hot plasma in the GC region (Koyama et 
al. 1996). Although the X-ray continuum could again be ex- 
plained by thermal free-free emission at kT > 10 keV, the 
simultaneous appearance in the spectra of Ka lines from 
helium-like and hydrogen-like ions of various elements sug- 
gests that the very hot plasma is either multi-temperature 
or out of ionization equilibrium. Its electron density was 



estimated by Koyama et al. (1996) at ^ (0.3 — 0.4) cm"-^. 
When integrated over an ellipsoid Gaussian distribution 
with FWHM size ~ (270 pc)^ x 150 pc (from Yamauchi 
et al. 1990), this electron density implies a total H+ mass 
~ (8 - 11) X 10^ Mq, 58% of which hes inside r = 150 pc. 
From this, we may conclude that the very hot gas encloses 
~ 7% - 10% of the ionized gas, and hence ~ 0.3% - 0.5% 
of all the interstellar gas, inside r = 150 pc. 



3. Dynamical models 

The observed non-axisymmetric distribution and non- 
circular motions of interstellar gas near the GC can most 
likely be attributed to the presence of a Galactic bar, and to 
a large extent, they can be explained in terms of the theo- 
retical properties of stable closed orbits in the gravitational 
potential of a barred galaxy (see Binney & Merrifield 1998 
for a brief review). Both analytical and numerical calcula- 
tions of particle orbits in a barred gravitational potential 
have brought to light the existence of two families of stable 
closed orbits. Although the exact status of these two fami- 
lies depends on the shape and strength of the bar and can 
become rather complex, for our purposes it suffices to be 
aware of the basic typical situation. Typically, the orbits 
prevailing inside the bar's inner Lindblad resonance (ILR) 
are the so-called X2 orbits, elongated perpendicular to the 
bar's major axis, and the orbits prevailing outside the ILR 
are the xi orbits, generally elongated along the bar - except 
between the bar's corotation radius and its outer Lindblad 
resonance (OLR), where most closed orbits are unstable 
(see, e.g., Contopoulos & Mertzanides 19770 Contopoulos 
& Papayannopoulos 1980; Athanassoula 1992a). The cUip- 
ticity of both families of orbits increases toward the ILR, 
and on the whole xi orbits are more elongated than X2 
orbits. 

Binney et al. (1991) proposed a dynamical model of in- 
terstellar gas in the innermost (2 — 3) kpc of the Galaxy 
meant to reproduce at best the observed morphology and 
kinematics of CO, CS and Hi emissions. The basic scenario 
underlying their model is the following: Interstellar gas sub- 
ject to the gravitational potential of the bar tends to settle 
onto stable closed orbits, although due to diverse dissipa- 
tion processes, it gradually drifts inward through a sequence 
of decreasing-energy orbits. Gas starting off outside the ILR 
travels along xi orbits until its inward drift brings it to a 
cusped xi orbit, interior to which xi orbits become self- 
intersecting at their apocenters. Shocks and material col- 
lisions along the cusped and outermost self-intersecting xi 
orbits cause the gas to abruptly lose energy and angular 
momentum and to drop onto the lower-energy X2 orbits, 
along which it continues its spiraling motion toward the 
GC. Since X2 orbits are overall rounder than xi orbits, gas 
loses energy and drifts inward less rapidly in the "x2 disk" , 
which therefore attains higher densities. 

Binney et al. (1991) interpreted the observational 
longitude- velocity {I, v) diagrams of CO and Hi emissions in 
terms of the theoretical (/, v) curves generated by non self- 
intersecting xi orbits and by X2 orbits. The cornerstone of 
their discussion is the parallelogram that encloses most of 
the emission in the CO (l,v) diagram (this parallelogram 



It is interesting to note that the xi versus xi terminology in 
use today is opposite to that initially introduced by Contopoulos 
& Mertzanides (1977) in the case of a weak bar. 
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corresponds to the 180-pc molecular ring discussed earlier, 
which, as we already mentioned, is also known as the "ex- 
panding molecular ring"). Binney et al. argued that the 
parallelogram could be identified with the (l.v) trace of 
the cuspcd xi orbit, thereby placing the GC molecular disk 
(i.e., the GC disk population of dense molecular clouds) on 
X2 orbits. No such parallelogram appears in the Hi {l,v) 
diagram, whose structure suggests instead that atomic gas 
extends well into the region of non self-intersecting xi or- 
bits. These considerations led Binney et al. to propose that 
a substantial fraction of the interstellar gas turns molecular 
at the strong shocks present in the vicinity of the cusped 
Xi orbit. Farther out, the observed gap in the interstellar 
gas distribution between radii ^1.5 kpc and ^ 3 kpc could 
presumably be explained by the generally unstable char- 
acter of particle orbits in a region extending radially from 
somewhere inside corotation to the OLR and by the result- 
ing expulsion of gas away from the unstable region. The 
accumulation of expelled gas just outside the OLR would 
then lie at the root of the Galactic disk molecular ring. 

Detailed quantitative comparisons between the observa- 
tional {I, v) diagrams and the theoretical {I, v) curves en- 
abled Binney et al. (1991) to place constraints on the pa- 
rameters of their model. The best constrained parameter 
is the inclination angle of the bar's major axis to the line 
of sight, which was found to be ^bar — 16° (within a few 
degrees) toward positive longitudes. Much more uncertain 
are the bar's dimensions; Binney et al. arbitrarily adopted 
ao = 1.2 kpc for the bar's semi-major axis (more exactly, 
the radius along the major axis beyond which the matter's 
density falls off very steeply) and 4:3:3 for the axis ratios, 
but the bar could easily be longer and more elongated. The 
bar's corotation radius was arbitrarily set to rcR = 2.4 kpc, 
consistent with the notions that the bar terminates some- 
what inside corotation (Contopoulos & Papayannopoulos 
1980; Sellwood & Sparke 1988; Athanassoula 1992b) and 
that the region around corotation is depopulated of gas. 

For all its successes, Binney et al.'s (1991) model leaves 
a number of important observational features unexplained. 
Most notably, the longitudinal lopsidedness in the observed 
CO distribution is much more pronounced than the asym- 
metry expected from the bar's inclination toward positive 
longitudes. In addition, Binney et al. noted a significant dis- 
crepancy in the velocity scales of the observational CO and 
Hi (Z, v) diagrams, when interpreted in the framework of 
their model. They attributed this discrepancy to the neglect 
of gas-dynamical effects in their particle-orbit approach, ar- 
guing that hydrodynamical forces at the strong shocks re- 
sponsible for the transition from Xi to X2 orbits near the 
ILR lead to large departures from ballistic trajectories. 

This argument prompted Jenkins & Binney (1994) to 
simulate interstellar gas flows in Binney et al.'s (1991) bar 
model. They tested two radically different two-dimensional, 
sticky-particle schemes, which neglect self-gravity, but al- 
low for both accretion of matter shed by bulge stars and 
removal of clouds from regions above a given density thresh- 
old. Both schemes yield similar results, which qualitatively 
confirm the dynamical scenario envisioned by Binney et 
al., with one important difference. Here it is found that gas 
parcels approaching (but still outside) the cusped xi or- 
bit tumble onto X2 orbits located deeper inside, so that a 
nearly empty zone develops near the cusped orbit between 
the "xi ring" and the "2:2 disk" . Only an external source of 
matter, such as gas loss by bulge stars, can keep this tran- 



sition zone appreciably populated. For future reference, the 
ILR lies at tilr — 350 pc, the cusped orbit has a half-size 
~ 650 pc X 150 pc, and for reasonable values of the numer- 
ical parameters, the inner boundary of the xi ring has a 
half-size ~ (750 ± 50) pc x (400 ± 50) pc and the X2 disk a 
half-size ~ (280 ± 20) pc x (180 ± 20) pc. 

Regrettably, Jenkins & Binney's (1994) simulations fail 
to resolve the main problems with Binney et al.'s (1991) 
model, namely, its inability to reproduce the pronounced 
lopsidedness in the observed CO distribution and its in- 
ability to simultaneously explain the observational CO and 
Hi {l,v) diagrams. Jenkins & Binney suggested that part 
of the blame could probably be laid on the limitations of 
their rnimerical schemes, which are restricted to two di- 
mensions and omit self-gravity. They also pointed out that 
the unsteadiness of gas flows near the ILR could be partly 
(though certainly not entirely) responsible for the observed 
asymmetry in CO. Finally, they raised the possibility that 
CO and Hi line emissions might not be reliable tracers of 
interstellar gas mass. 

Subsequently, many hydrodynamical simulations of gas 
flows in barred gravitational potentials were performed 
with various numerical techniques. Englmaier & Gerhard 
(1999) employed a two-dimensional smoothed-particle hy- 
drodynamics code to follow interstellar gas flows in the 
Galactic gravitational potential deduced from the depro- 
jectcd and dust-corrected COBE/DIRBE near-infrared lu- 
minosity distribution (see Binney et al. 1997), under the as- 
sumption of a constant mass-to-light ratio. The COBE data 
imply that the bar's inclination angle is restricted to the in- 
terval 15° - 35°, with a preferred value 0bar = 20°. With 
this value, the bar has a semi-major axis a ~ 2 kpc and axis 
ratios ~ 5 : 3 : 2, it is surrounded by a thin elliptical disk 
with radius 3.5 kpc along the major axis and ~ 2 kpc 
along the minor axis, and corotation occurs at rcR, ~ 3 kpc 
(assuming tq = 8 kpc; Binney et al. 1997). Englmaier & 
Gerhard (1999) supplemented the gravitational potential 
deduced from the COBE near-infrared luminosity with the 
contributions from a central cusp and (for some runs) an 
outer dark halo. They presented their simulation results 
in the form of face-on surface-density maps of interstellar 
gas in the inner Galaxy (r < r©) together with the associ- 
ated (Z, v) diagrams. The maps clearly feature a four-armed 
spiral pattern outside corotation and two pairs of arms in- 
side corotation, each pair emanating approximately from 
one end of the bar. Interior to the two innermost arms, the 
vicinity of the cusped xi orbit becomes rapidly depleted, as 
gas parcels reaching the cusped orbit's shock quickly fall in 
along the shock ridges onto the central X2 disk. The latter 
stands out in the surface-density maps as a compact disk 
of radius ~ 150 pc. 

Bissantz et al. (2003) ran similar hydrodynamical sim- 
ulations based on the improved near-infrared luminosity 
distribution models of Bissantz & Gerhard (2002). In the 
best-flt model that they adopted, the bar's inclination an- 
gle is again ^bar = 20°, but its semi-major axis is now 
a ~ 1.75 kpc, its axis ratios ~ 5 : 2 : 1.5 and its corota- 
tion radius rcn — 3.4 kpc (assuming again = 8 kpc). 
Bissantz et al. (2003) examined the gas dynamics in the 
bar region more closely. As a starting point, they computed 
closed xi and X2 orbits in their model gravitational poten- 
tial and found that the outermost X2 orbit does not reach 
out to the cusped xi orbit. This, they argued, could be the 
reason why a gap opens up in the simulations between the 



12 



K. Ferriere et al.: Spatial distribution of interstellar gas in the Galactic bulge 



cusped orbit and the X2 disk. According to their calcula- 
tions, the cusped orbit has a half-size ~ 1.35 kpc x 280 pc 
(approximately twice as large as in Jenkins & Binney 
1994) and the outermost closed X2 orbit has a half-size 
~ 175 pc X 85 pc (nearly twice smaller than the half-size of 
the X2 disk in Jenkins & Binney 1994). 

All the above simulations have two important limita- 
tions: they are two-dimensional and, more critically, the 
Galactic gravitational potential is prescribed at the out- 
set and taken to be time-independent and point-symmetric 
with respect to the GC. A major improvement on both 
counts was brought by Fux (1999), who performed three- 
dimensional, self-consistent composite simulations of the in- 
ner Galaxy, in which the stellar bar and the interstellar gas 
are evolved in concert, with iV-body and smoothed-particle 
hydrodynamics codes, respectively. All gravitational inter- 
actions (including self-gravity) are taken into account, and 
no stationarity or symmetry condition is imposed. The sim- 
ulated face-on surface-density maps are found to fluctuate 
in time. The stellar bar varies in shape and its center of 
mass oscillates around the GC with an amplitude of several 
100 pc. The interstellar gas flow is unsteady and asymmet- 
ric; it develops transient spiral arms and shock fronts as 
well as a ring or disk of X2 orbits which closely follows the 
bar's center of mass. 

Fux (1999) constrained the location of the Sun relative 
to the bar with the help of dust-corrected COBE/DIRBE 
K-band data. Moreover, he constructed movies of synthetic 
(^, v) diagrams, compared them to the observational CO 
and Hi {I, v) diagrams and selected the best-matching snap- 
shots. He then used the corresponding gas flow models to 
interpret the main features of the observational (/, v) dia- 
grams. The two preferred models have 6'bar — 25° (com- 
patible with the COBE data) and rcn — 4.5 kpc. In the 
model that Fux singled out for attention (the one that gives 
the best overall agreement with the {l,v) observations), 
the cusped xi orbit has a half-size ~ 3.2 kpc x 1.3 kpc, 
which is considerably larger than in the simulations de- 
scribed above. With such large dimensions, the cusped orbit 
cannot possibly correspond to the parallelogram in the ob- 
servational CO {I, v) diagram, as suggested by Binney et al. 
(1991). Instead, Fux's conjecture is that the parallelogram 
could represent interstellar gas streams originating from the 
cusped orbit or from dustlanes (associated with "off-axis 
shocks") connecting the cusped orbit to the X2 ring/disk, 
and grazing the X2 ring/disk. Aside from this new inter- 
pretation, Fux pointed out that the bar parameters cannot 
be estimated with any confidence from the best-matching 
models, and he devised a model-independent method rely- 
ing on purely geometric considerations. This method yields 
^bar = 25° ± 4° for the bar's inclination angle, a ~ 3.2 kpc 
for its semi-major axis, a : & ~ 5 : 3 for its face-on axis 
ratio and rcn — (4.0 ± 0.5) kpc for its corotation radius 
(assuming = 8 kpc). 

A great merit of Fux's (1999) study is that it under- 
scores the inherent non-stationarity of the Galactic bar, 
with its large- amplitude oscillations around the GC, and 
the resulting unsteadiness and off-centering of the inter- 
stellar gas distribution. In this manner, it offers a natural 
explanation for the observed longitudinal asymmetry in the 
CO intensity maps. Furthermore, since the cusped xi orbit 
is no longer associated with the parallelogram in the CO 
(/, v) diagram, the difference in the velocity scales of the 
observational CO and Hi (/, v) diagrams no longer poses a 



problem. On the other hand, Fux's simulations raise their 
own set of questions, as the large values obtained for the 
corotation radius and for the dimensions of the cusped or- 
bit seem difficult to reconcile with the observed gas deficit 
between r ~ 1.5 kpc and 3 kpc and with the Galactic disk 
molecular ring starting at r ~ 3.5 kpc. 

Another, more subtle aspect of the interstellar gas mor- 
phology near the GC concerns the orientation of the cen- 
tral X2 ring/disk with respect to the bar's major axis. 
We know that closed particle (i.e., coUisionless) orbits in 
a barred gravitational potential are elongated either par- 
allel or perpendicular to the bar's major axis. Typically, 
they are elongated parallel to the bar between the ILR and 
corotation and outside the OLR, and perpendicular to the 
bar inside the ILR (where the central X2 ring/disk resides) 
and between corotation and the OLR (e.g., Contopoulos & 
Mertzanides 1977; Contopoulos & Papayannopoulos 1980). 
This is not true of closed gas (i.e., coUisional) orbits. Once 
pressure and viscous forces are taken into account, the 
90° jumps in orbit orientation at the bar's resonances are 
completely smoothed out. Wada (1994), who studied the 
case of a weak bar with an analytical damped-oscillator 
model, found that the orbit inclination angle to the bar's 
major axis, A9, varies gradually with radius, such that 
0° < A9 < 45° between the ILR and corotation and outside 
the OLR, and 45° < A0 < 90° inside the ILR and between 
corotation and the OLR. Hence, X2 orbits lead the bar by 
A6 > 45°, and the whole X2 ring/disk should appear in- 
clined to the bar by a little more than 45°. This is exactly 
what a close scrutiny of the face-on maps of Englmaier & 
Gerhard (1999) and Fux (1999) shows. In addition, if the 
bar is itself inclined by 6'bar — 20° — 25° to the line of 
sight (Englmaier & Gerhard 1999; Fux 1999; Bissantz et al. 
2003), the X2 ring/disk should be inclined by ~ 70° to the 
line of sight, in excellent agreement with the observational 
conclusions of Sawada et al. (2004). Let us also note that 
the gradual variation in orbit inclination with radius au- 
tomatically gives rise to spiral density enhancements. One 
could surmise that these enhancements are related to the 
molecular arms discussed by Sofue (1995a) (see section l^^ . 

4. Our model for the gas distribution 

It is clear that our observational knowledge of the interstel- 
lar gas distribution in the GB is rather patchy. Although 
theoretical models can probably bridge some of the gaps, 
there remain important murky areas. In this section, we 
gather the relevant observational results (see section[2|) that 
we deem the most trustworthy, we complement them with 
some theoretical predictions from gas dynamical models 
(see section [3]) and we try to piece everything together into 
a coherent picture of the interstellar gas distribution in the 
innermost 3 kpc of our Galaxy. 

4.1. The central molecular zone 

There have been several attempts to deproject sky maps of 
molecular line emission in the CMZ. In most cases, a gas 
kinematic model is appealed to in order to transform the 
measured line-of-sight velocity into line-of-sight distance. 
Unfortunately, our poor knowledge of the true gas kine- 
matics near the GC renders this kind of method unreliable. 

A notable exception is the work of Sawada et al. (2004), 
who derived a face-on map of the molecular gas without 
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any kinematic assumption, on the sole basis of observa- 
tional (emission versus absorption) data (see section \2.2^ . 
For this reason, their face-on map will serve as the first 
building block of our model. Thus, we assume that the 
CMZ projected onto the Galactic plane has the shape 
of a 500 pc X 200 pc ellipse inclined by 70° to the line 
of sight toward positive longitudes. From Sawada et al.'s 
Figure 10a, we estimate that the ellipse is centered on 
(xcUc) — (—50 pc, 50 pc)0 The peak in surface den- 
sity appears displaced eastward from the ellipse center, to 
~ (—30 pc, 100 pc), but in view of the large uncertainties 
involved in the deprojection, the apparent peak position 
should not be taken at face value. It is even possible that 
the molecular gas is in fact distributed in a ring rather 
than a disk (e.g., Sofue 1995a). While we do not rule out 
this possibility, we find the present observational evidence 
inconclusive. Theoretical studies are not of great help ei- 
ther. Accepting the idea that the CMZ corresponds to the 
domain of X2 orbits (see section [3]) , we note that both X2 
disks and X2 rings have come up in hydrodynamical simu- 
lations (e.g., Jenkins & Binney 1994; Englmaier & Gerhard 
1999; Fux 1999). 

Here, as a simple compromise, we adopt a nearly flat 
(i.e., between peaked and ring- like) horizontal distribution 
of the type utilized by Launhardt et al. (2002) for their 
warm inner disk (Equation [T]) . For convenience, we intro- 
duce horizontal coordinates in the CMZ frame, {X,Y), de- 
fined such that X is along the major axis and Y along the 
minor axis. The CMZ coordinates {X, Y) are then related 
to the Galactic coordinates {x, y) through 



X = [x - Xc) cos 9c + [y -Vc) sin 9c 
Y = -{x - Xc) sm9c + (y - yc) cost 



(15) 



with Xc = —50 pc, yc = 50 pc and 9c = 70° (see Figured]). 
Adjusting the horizontal profile from Equation [1] to an el- 
liptical disk with semi-major axis 250 pc and axis ratio 2.5 
gives for the H2 space-averaged density in the CMZ 



("-H2/CMZ oc cxp 



V^2 + (2.5y)2-Xc 



, (16) 



where Xmax = 250 pc, Xc = Xniax/2 = 125 pc and Lc = 
^max/(2 (In 2)1/4) ^ Evidently, Equation [H] can 

easily be adapted to a larger/smaller CMZ or to a more/less 
elongated CMZ by increasing/decreasing the value of X^ax 
or the prefactor of F, respectively. It can also be adapted to 
a ring morphology by increasing Xc/Xj^ax and decreasing 
accordingly. 

Sawada et al. (2004) provided only a face-on view of 
the molecular gas, with no information on its vertical dis- 
tribution or density scale. Regarding the vertical distribu- 
tion, one has to distinguish between the GC molecular disk 
and the surrounding 180-pc molecular ring (see section l^?^ . 
The GC molecular disk appears planar and closely aligned 
with the Galactic plane (Heiligman 1987; Bally et al. 1988; 
Burton & Liszt 1992; Oka et al. 1998) . Its FWHM thick- 
ness is typically ~ 15 — 30 pc, with a tendency to increase 
outward (up to 50 pc) as well as toward massive molecu- 
lar complexes (up to '--^ 30 — 60 pc) (Heiligman 1987; Bally 



^ The {x, y, z) coordinate system used by Sawada et al. (2004) 
is equivalent to our [x, y, z) system (see Figure [ij rotated clock- 
wise by 90°. 




Fig. 2. Face-on view of our CMZ. {x,y,z) are the 
Galactocentric coordinates and {X, Y) the CMZ horizon- 
tal coordinates. The CMZ is elliptical, 500 pc x 200 pc in 
size, centered on {xc,yc) = (—50 pc, 50 pc) and inclined by 
9c = 70° to the hue of sight. 



et al. 1988; Sofue 1995a; Burton & Liszt 1992). Here, follow- 
ing Burton & Liszt (1992), we ignore the spatial variations 
of the disk thickness, which we set to 30 pc, and we assume 
a Gaussian vertical distribution: 



(nH2)cMZ oc exp 



(17) 



with He = (30 pc)/(2 \/W2) ~ 18 pc. The 180-pc molecu- 
lar ring is thicker than the GC molecular disk (Bally et al. 
1988; Sofue 1995b), consistent with the outward increase 
in the GC molecular disk's thickness noted above, with the 
flaring of the H2 disk of Burton & Liszt (1992) and with 
the greater scale height obtained by Sanders et al. (1984) 
for the H2 GB disk as a whole (r < 1.3 kpc). In addi- 
tion, the 180-pc molecular ring is tilted out of the Galactic 
plane, counterclockwise by ~ 5° — 9° (Heiligman 1987). It 
is also inclined to the line of sight, as implied by the fact 
that its negative-velocity portion (near side) lies at nega- 
tive latitudes, while its positive-velocity portion (far side) 
lies at positive latitudes (Heiligman 1987; Bally et al. 1988). 
Heiligman (1987) estimated an inclination angle ~ 85° to 
the plane of the sky, whereas from Figure 11 of Bally et al. 
(1988) together with the attendant discussion we estimate 
an inclination angle ~ 80°. It is noteworthy that both the 
tilt and the inclination of the 180-pc molecular ring are in 
the same sense as those of the GB disk (see section 12.21 for 
the molecular GB disk and section [^31 for the Hi GB disk). 
In any event, since the 180-pc molecular ring actually con- 
tains but a small fraction of the CMZ mass, we will not 
treat it separately from the GC molecular disk. 

We now turn to the H2 density scale, which can be in- 
ferred from the total H2 mass of the CMZ. Our H2 mass 
estimates (assuming Xqq = 5 x 10^^ cm~^ K"^ km"^ s) 
for various regions of the GB are summarized in Table [4l 
Clearly, none of the listed region corresponds exactly to our 
CMZ, whose projected half-size is 250 pc x sin 70° ~ 235 pc. 
The closest is the region r < 200 pc studied by Sofue 
(1995ab), which contains ^ 1.3 x lO'^ Mq of H2 in the 
GC molecular disk plus - 2.7 x 10^ Mq of H2 in the 180- 
pc molecular ring, i.e., a total H2 mass ^ 1.6 x 10^ Mq. 
From this, we estimate that our CMZ has an H2 mass 
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^ 1.9 X 10^ Mq. This value is in good agreement with 
the H2 mass of the GB region < 375 pc being ^ 
(4.2 - 4.4) X 10^ Mq (Oka et al. 1998; Sodroski et al. 1995; 
see Table U) and with the total hydrogen mass of the region 
r < 230 pc being - 2 x lO'^ Mq (Launhardt et al. 2002; see 
Table [1]). 

The H2 density scale can now be obtained by inte- 
grating the density model described by Equations \W\ and 
1171 over the relevant region and equating the integral to 
our estimated H2 mass of 1.9 x 10'' Mq. The result is 
(('^H2)cMz)max — 150 cm~'^, SO that the H2 space-averaged 
density in the CMZ can be written as 



space-averaged density in the CMZ is given by 



(»^H2)cMZ = (150 cm ^) exp 



VX2 + (2.5r)2-Xe 



X exp 



(18) 



where = 125 pc, = 137 pc, = 18 pc and {X,Y) 
are the CMZ horizontal coordinates defined by EguationfTSl 
On the plane of the sky, the CMZ appears as a 474 pc x 30 pc 
(size at half-maximum density) ellipse, displaced eastward 
by 50 pc, so that it extends out to r± = 287 pc at positive 
longitudes and rj_ = 187 pc at negative longitudes (see 
Figure H^) . 

For the atomic gas, all the face-on maps deduced from 
21-cm observations rely on a number of kinematic assump- 
tions, which render them questionable. Therefore, we prefer 
to take an alternative approach, based on our more solid 
knowledge of the molecular gas distribution and on the 
premise that atomic and molecular gases are similarly dis- 
tributed in space (as argued by Liszt & Burton 1996). The 
only obvious difference between both gases concerns the 
thickness of their layer, which is ~ 30 pc for H2 (see above 
Equation [17]) and ~ 90 pc for Hi (Rohlfs & Braunsfurth 
1982; see section [^75)1 . Thus, we adopt for Hi the same hor- 
izontal and vertical profiles as for H2 (Equations [11] and 
[T7[ respectively), with this difference that He ^ H'^ — 
(90 pc)/(2%/h^) ~ 54 pc. 

The Hi mass of the CMZ is very uncertain. As we men- 
tioned in section [2T3I current Hi mass estimates differ by up 
to a factor ~ 5. Here, we favor the mass estimates of Rohlfs 
& Braunsfurth (1982), who directly focused on a region 
comparable in size to our CMZ - as opposed to Burton & 
Liszt (1978, 1993) and Liszt & Burton (1980), who studied 
the GB disk as a whole. Rohlfs & Braunsfurth (1982) de- 
rived an Hi mass ~ 3.1 x 10® Mq for the region r < 300 pc, 
which translates into an Hi mass ~ 1.4 x 10® Mq for the 
region r < 200 pc. Since the H2 mass in the latter region 
is ~ 1.6 X 10'' Mq (Sofue 1995ab; see paragraph preceding 
Equation [TS]), we may conclude that the Hi mass amounts 
to 8.8% of the H2 mass. 

With three times the scale height and 8.8% of the mass 
of the molecular gas, the atomic gas reaches only 2.9% of its 
maximum mass density, or 5.8% of its maximum number 



\"Hi/CMZ 



cm ) exp 



X exp 



(19) 



with HI, — 54 pc and Xc, Lc, {X,Y) as defined below 
Equation [TBI fsee Figure [1|d). 

4.2. The "holed" GB disk 



density, i.e., ((nHi)cMz) 



8.8 cm""*. Altogether, the Hi 



We saw in section 14.11 that reliable observational informa- 
tion on the interstellar gas distribution inside the CMZ is 
scanty. It turns out that the observational situation is even 
worse outside the CMZ. It is true that a number of tilted 
GB disk models have been constructed on the basis of CO 
and 21-cm observations (see sections 12.21 and 12.31) . but all 
of them involve unverified assumptions, both on the overall 
geometry (e.g., axial symmetry) and on the gas kinematics. 

The model that we consider the most plausible, despite 
its inherent uncertainties, is that of Liszt & Burton (1980), 
because it is the only model with an elongated disk, which 
makes it the easiest to reconcile with our theoretical under- 
standing of gas dynamics in the GB (see section[3]). Indeed, 
the GB disk outside the CMZ is believed to correspond to 
the domain of non self-intersecting xi orbits, which form 
an elongated, radially thick ring around the X2 disk/ring 
(itself identified with the CMZ). In some hydrodynamical 
simulations, the inner boundary of the xi ring grazes the 
X2 disk/ring, while in others, a gap opens up between both. 

Here, we start from Liszt & Burton's (1980) GB disk, 
namely, an elliptical disk with semi-major axis 1.6 kpc and 
axis ratio 3.1, and we introduce a hole in the middle just 
large enough to enclose the CMZ. It can be checked that a 
hole with the same shape and half the size of the entire GB 
disk is sufficient for our purpose. Then the inner boundary 
of the "holed" GB disk has a half-size of 800 pc x 258 pc, 
in reasonably good agreement with the results of hydrody- 
namical simulations. For instance, Jenkins & Binney (1994) 
found that the inner boundary of the xi ring has a half-size 
~ (750±50) PCX (400±50) pc. The xi ring is more difficult 
to make out in the larger-scale simulated maps of Englmaier 
& Gerhard (1999) and of Bissantz et al. (2003), but its in- 
ner boundary clearly appears more elongated, with an axis 
ratio possibly close to our adopted value of 3.1. Let us also 
point out that our GB disk is more elongated than our CMZ 
(whose axis ratio is 2.5; see section Wl} . This is perfectly 
consistent not only with gas fiow simulations, but also with 
particle orbit calculations, which predict a greater elliptic- 
ity for xi orbits than for X2 orbits (see first paragraph of 
section [3]) . 

According to Liszt & Burton (1980), the GB disk is 
tilted by a = 13?5 out of the Galactic plane [counterclock- 
wise rotation by a about the x-axis, whereby (x, y, z) 
{x,y', z')] and inchned by j = 70° to the plane of the sky 
["front-downward" rotation by f3 = 90° — i about the y'- 
axis, whereby {x,y\z') {x" ,y' , z")], and the GB disk's 
major axis forms an angle = 48?5 to the a;"-axis [mea- 
sured clockwise about the 2;"-axis] . If we denote the spatial 
coordinates in the GB disk frame by {X,y,Z), with X 
along the major axis, y along the minor axis and Z along 
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the northern normal, the transformation from Galactic co- 
ordinates to GB disk coordinates reads 



r X 



y 



z = 



X cos/9 cos^d 

— y (sin a sin/3 cos6'd — cosa sin^d) 

— z (cos a sin/? cos6'd + sinQ; sin6'd) 
—X cos/3 sin^d 

+ y (sin a sin/3 sin^d + coso; cos^d) 
+ z (cos a sin/3 sin 0d^ sin a cos^d) 
X sin/3 

+ y sin a cos/3 
+ z cos a cos/3 , 



(20) 



with a = 13?5, /3 = 20° and ^d = 48?5 (see Figure IS]). 

Now that we have specified the horizontal shape and 
dimensions as well as the orientation of our holed GB disk, 
we may fit to it horizontal and vertical profiles of the same 
type as for our CMZ (Equations [TBI and [Tfl respectively). 
Again, we assume that molecular and atomic gases are sim- 
ilarly distributed in space, except for a possible difference 
in scale height. Under these conditions, we can write 



(?iH2/disk oc exp 



l»^Hi)disk oc exp 



Ld 



with A'max = 1-6 kpC, A'miii = 0.8 kpc, Xd = (A'max 

A'„i„)/2 - 1.2 kpc and = {X,,,,^ - X^,^)l{2 (ln2)i/4) 
438 pc, and 



(»^H2/disk oc exp 



^iHi/disk oc exp 



Z 



H'd 



For the scale height, we adopt Hd = (70 pc)/(2vln2) ~ 
42 pc for H2 (Sanders et al. 1984) and H'^ = V2 (85 pc) ~ 
120 pc for Hi (Liszt & Burton 1980). Like in the CMZ, the 
atomic layer is about three times thicker than the molecular 
layer. Moreover, each layer is ~ 2.2 — 2.3 times thicker in 
the holed GB disk than in the CMZ. 

Our density model can again be normalized with the 
help of the total H2 and Hi masses in the considered re- 
gion. The entire GB disk contains ~ 5.3 x 10^ Mq of H2 
(Sanders et al. 1984; see Table S]) and - 5.2 x 10^ Mq of 
Hi (Liszt & Burton 1980; see Table[5]), while our CMZ con- 
tains - 1.9 X 10^ Mq of H2 and - 1.7 x 10*^ Mq of Hi 
(see section BTTj) . By subtraction, we find that our holed 
GB disk has an H2 mass ~ 3.4 x 10'' Mq and an Hi 
mass ~ 3.5 x 10^ Mq. We are fully aware of the pitfalls 
of this mass estimation, especially as the masses involved 
in the subtraction procedure have considerable uncertainty. 
Particularly uncertain is the H2 mass of the entire GB disk. 
The only estimate available to us is the value from Sanders 
et al. (1984). But this value applies to an axisymmetric disk 
of radius 1.3 kpc having twice the surface area of our el- 
liptical GB disk, which means that it could overestimate 
the true mass by up to a factor ~ 2. On the other hand, 
the comparatively large values obtained by Sodroski et al. 




Fig. 3. Three-dimensional view of our holed GB disk. 
{x,y,z) are the Galactocentric coordinates, {X,y,Z) the 
GB disk coordinates, and {x,y',z'), {x",y',z") the auxil- 
iary coordinates defined above Equation [20l The holed GB 
disk is elliptical, 3.2 kpc long and (3.2 kpc)/3.1 ~ 1.03 kpc 
wide, and its inner boundary has half the dimensions of its 
outer boundary. Relative to the Galactic plane, its principal 
plane is tilted by a = 13? 5 about the x-axis and inclined 
by /3 = 20° about the y'-axis, and its major axis is oriented 
at 6d = 48? 5 to the x"-axis. 



(1995) and by Oka et al. (1998) for a smaller region (see 
Table |4]) suggest instead that Sanders et al.'s (1984) value 
could underestimate the true mass by roughly the same fac- 
tor. Hence, the two sources of error tend to counteract each 
other. In addition, two facts allow us to place some confi- 
dence in our results. First, the H2 and Hi masses of our 
holed GB disk are slightly smaller than what they would 
be if the gas were uniformly distributed throughout the GB 
disk (- 4.0 X 10"^ Mq and - 3.9 x 10^ Mq, respectively). 
Second, Hi represents 8.8% by mass of H2 in the CMZ and 
10.3% in the holed GB disk. This statement supports Liszt 
& Burton's (1996) claim that both gases have similar spa- 
tial distributions, not only in the CMZ and in the holed 
GB disk separately (as we already assumed earlier), but 
also throughout the entire GB disk. It is also reassuring to 
obtain a slightly lower Hi mass fraction in the CMZ. 

If we now equate the space integral of Eauation l21l times 
Equation [22] to our estimated H2 and Hi masses, we find 
((nHa) disk) max ^ 4.8 cm.-^ and ((nHi)disk)max 0.34 cm-3. 
Hence, the holed GB disk has an H2 space-averaged density 



(?^H2)disk = (4.8 cm ) exp 



y/X^ + i3.iy)^-Xd^ 

Ld I 



X exp 



Z 
Ifd 



(23) 



and an Hi space-averaged density 



(nHi)disk (0.34 cm ^) exp 



' ^ X^ + {2,.iyY - Xd' 
Ld 



X exp 



H'd 



(24) 
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Fig. 4. Projection of the CMZ (bright area) and the holed 
GB disk (fainter area) onto the plane of the sky: (a) molec- 
ular gas (see Equations [18] and [23|) ; (b) atomic gas (see 
Equations fT9l and [24)l . The apparent sizes are a little larger 
than the sizes at half-maximum density, because of the log- 
arithmic scale used in the projection. In contrast to the 
CMZ, which is truly displaced to the left, the GB disk is 
symmetric with respect to the GC, and the only reason 
why it appears more extended on the left side is because 
its positive-longitude portion lies closer to us. 



with Xd = 1.2 kpc, Ld = 438 pc, = 42 pc and H'^ = 
120 pc. On the plane of the sky, the GB disk extends out to 
r_L = 1.14 kpc (radius at half-maximum density) on each 
side of the GC (see Figure S]). Projected onto the Galactic 
plane, it has the shape of a 2.94 kpc x 1.02 kpc (FWHM 
size) ellipse inclined clockwise by 47? 6 to the line of sight 
(see Figure [5]). This inclination angle is greater than that 
typically found for the Galactic stellar bar (6'bar — 15°— 35°; 
see section [3]) , but it is in good agreement with the value 
^bar = 44° ±10° recently obtained by Benjamin et al. (2005) 
from the GLIMPSE Point Source Catalog. 



4.3. The ionized component 

The best available model for the spatial distribution of in- 
terstellar free electrons in the GB is the NE2001 model of 
Cordes & Lazio (2002) presented in section [2T4l According 
to this model, the total mass of interstellar ionized hydro- 
gen in the region r < 3 kpc is (7.3x lO'^ MQ)/(l-h0.2 /him), 




y [kpc] 

Fig. 5. Projection of the CMZ (bright area) and the holed 
GB disk (fainter area) onto the Galactic plane. Displayed 
here is the H2 map (from Equations [18] and [23]) . The Hi 
map (from Equations [T9l and [24| looks identical, except for 
this hardly noticeable difference that the GB-disk-to-CMZ 
luminosity ratio is slightly greater. For the same reason as 
in Figure H] the apparent sizes are a little larger than the 
sizes at half-maximum density. 



where /him is the fraction of ionized gas belonging to the 
hot medium (see Table HJ). The mass of hot 11+ in the same 
region can be estimated from Almy et al.'s (2000) model 
(neglecting the contribution from very hot 11+) at 1.2 x 
10"^ Mq (see Table [5]). It then follows that /him = 17% (or, 
equivalently, /wim = 83%) and that the total mass of 11+ 
inside 3 kpc is 7.1 x lO'^ Mq, divided between 5.9 x lO"^ M© 
in the WIM and 1.2 x lO"^ Mq in the HIM. Furthermore, 
from Equation [Tl] with /him = 17%, we gather that the 
H+ space-averaged density is given by (nH+) — 0.97 {ric). 
The partial contributions from the warm and hot ionized 
media are globally given by (nH+)wiM = /wim {nn+) and 
("-h+)him = /him (?iH+), respectively. For the WIM, which 
contributes a large 83% of the total H+ mass, we may rea- 
sonably assume that the above global relation remains ap- 
proximately valid locally. Owing to the large uncertainties 
in the exact spatial dependence of the density distributions, 
we feel that taking (»^h+)wim = /wim (f^H+) at all r is 
safer than subtracting (?ih+)him (which can be estimated 
independently; see next paragraph) from (nH+). In that 
case, the H+ space-averaged density of the WIM is simply 
("h+)wim — 0-80 ("c) or, in view of Equations 171 - [TUl 



("H+)wiM = (8-0 cm 3) 

x'^ + {y- y^f 



x < exp 



^3 



+ 0.009 exp 
+ 0.005 



r - L2 



exp 
2 



Hi 



sech 



2Li 



u{Li - r) 



Ho 



sech^ ( — 



Hi 
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with 2/3 = —10 pc, Z3 = —20 pc, L3 = 145 pc, H3 = 26 pc, 
L2 = 3.7 kpc, H2 = 140 pc, Li = 17 kpc and Hi = 950 pc. 

For the HIM, we turn again to Almy et al.'s 
(2000) hydrostatic-polytropic model. With (nH+)HiM — 
(ne)HiM/l-2, the H+ space-averaged density of the HIM di- 
rectly follows from Equation [HI 

(nH+)HiM = {(0.009 cm-3)2/3 _ 54 ^ iQ-i^ cm-" s^) 
X [<^(r,z)- 0(0,0)]}'"' , (26) 

where (j){r, z) is the Galactic gravitational potential of 
Wolfire et al. (1995) given by Equation [Ql 

For completeness, we also provide a description of the 
very hot medium detected by Koyama et al. (1989) and 
discussed at the end of section 12.41 Based on Yamauchi 
et al.'s (1990) and Koyama et al.'s (1996) studies, we as- 
sign to the VHIM an ellipsoid Gaussian density distribu- 
tion with central electron density 0.35 cm~^ and FWHM 
size (270 pc)^ x 150 pc. The major plane of the ellipsoid is 
seen edge-on from the Sun and tilted clockwise by 21° with 
respect to the Galactic plane, which leads us to introduce 
the coordinate transformation 



rj = y cos ttvh + z sin avh 
C — —y sinavh + z cos avh , 



(27) 



with Q!vh — 21°. The H+ space-averaged density of the 
VHIM, (nH+)vHiM = 
terms of the {rj, () coordinates as 

('^H+)vHiM = (0-29 cm~^) exp 

with Lyh = 162 pc and iJvh = 90 pc. 



c/vhim/1-2i can then be written in 



(28) 



4.4. Altogether 

The total space-averaged density of hydrogen nuclei in the 
interstellar GB is given by the sum of the partial contribu- 
tions from the molecular, atomic and ionized media: 



(n-n) = 2 (riHa 



where 



("-H2) = ("•H2)CMZ + ("'H2)disk 

(see Equations [TH] and [ 



("Hi) = ("-Hi)cMZ 

(see Equations [19] and [24]) and 



("H+>Hi 



\ "-Hi/disk 



(29) 

(30) 
(31) 
(32) 



(see Equations [23 [HI and [^5)1 . For an assumed total-to- 
hydrogen mass ratio of 1.453 (see beginning of section [2]), 
the total space-averaged mass density of interstellar gas in 
the interstellar GB is related to the above hydrogen space- 
averaged density through 



(p) - 1.453 TOP (nn) 



(33) 



where mp is the proton rest mass. 

Displayed in Figure [6| is the radial variation of 
the azimuthally-averaged column densities through the 




1.5 2 
r [kpc] 

Fig. 6. Azimuthally-averaged column densities through 
the Galactic disk of interstellar hydrogen nuclei in molec- 
ular form (solid line), atomic form (dot-dashed line) and 
ionized form (dotted line) and associated surface densities 
of total interstellar matter (assuming a total-to-hydrogen 
mass ratio of 1.453; see section [5]) as functions of Galactic 
radius. 



Galactic disk of the molecular, atomic and ionized gases. 
The vertical variation of their hydrogen space-averaged 
densities is shown in Figure [7| at three representative hor- 
izontal locations, namely, at the GC [{x,y) — (0,0)], at 
the center of the CMZ [{x,y) = (—50 pc, 50 pc)] and 
at the point of the GB disk (projected onto the Galactic 
plane) where its major axis intersects its "ridge" [r = 
1(1.47 kpc), = 47?6 (see below Equation [24]), i.e., 
(x, y) = (742 pc, 815 pc)]. The graphs speak for themselves. 
The molecular gas (solid line) is by far the most abundant. 
Its radial distribution exhibits a high peak at the origin cor- 
responding to the CMZ and a lower bump at r ~ 4.7 kpc 
corresponding to the holed GB disk. Along the vertical, it 
is confined to a thin layer, which is nearly centered on the 
Galactic midplane, except outside the CMZ, where it is cen- 
tered on the tilted midplane of the GB disk (see Figure [Tj:). 
The atomic gas (dot-dashed line) is about ten times less 
abundant than the molecular gas, it has virtually the same 
radial distribution and it occupies an approximately three 
times thicker layer. Finally, the ionized gas (dotted line) is 
the most rarefied, and it spreads much farther out in all 
directions. The slight rise in its radial distribution outside 
^ 1 kpc is due to the contribution from the Galactic disk. 



5. Conclusions 

In this paper, we took stock of the present observational 
status of interstellar gas in the Galactic bulge (defined as 
the region interior to r = 3 kpc) and we discussed in some 
detail existing theoretical models of interstellar gas dynam- 
ics near the GC. We then provided simple analytical ex- 
pressions for the hydrogen space-averaged densities of the 
different gas components that fit at best the observational 
data and are consistent with the theoretical predictions. 
All our expressions are written in parametric form, which 
will make them easy to update and adapt to future obser- 
vations showing other characteristics (shape, dimensions, 
orientation, mass) than those adopted here. Our results are 
summarized in section [4.41 
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Fig. 7. Space-averaged densities of interstellar hydrogen 
nuclei in molecular form (solid line), atomic form (dot- 
dashed line) and ionized form (dotted line) as functions 
of Galactic height: (a) at the GC; (b) at the center of the 
CMZ; (c) at the point of the GB disk (projected onto the 
Galactic plane) where its major axis intersects its "ridge" . 



We did not address the issue of where our modeled in- 
terstellar gas would appear in an (/,w) or full {l,b,v) plot. 
Doing so would require studying the gas dynamics in a 
specified non-axisymmetric Galactic gravitational potential 
with a given pattern speed. This is undoubtedly a step that 



should be taken in the future, especially as the (/, 6, v) space 
is the one that observers have direct access to. 

Aside from the space-averaged densities, we were able 
to gather some interesting observational information on the 
physical conditions (temperature and true density) in the 
molecular and ionized media. In section l^T^ we saw that the 
molecular medium probably contains a cold, high-density 
component with T ~ 50 K and ^ lO'^'^^"' cm~^, 
and a warm, low-density component with T ~ 150 K and 
riHa ~ 10^'^ cm~^. It is also possible that the temperature 
and true density in the molecular medium vary more-or-less 
continuously over broad ranges of values. In section we 
were led to describe the ionized medium as the superposi- 
tion of three components with radically different temper- 
atures: the WIM with T - 10"* K, the HIM with T ~ a 
few 10^ K and the VHIM with T > 10^ K. However, we 
were unable to obtain independent estimates for the asso- 
ciated true densities. As for the atomic medium, we may 
conjecture by analogy with the interstellar vicinity of the 
Sun that a cold and a warm phase co-exist, although we 
found little direct evidence to substantiate this hypothesis 
and even less direct quantitative information on the physi- 
cal conditions in these two putative phases. 

It would be highly desirable to estimate the filling fac- 
tors of the different media. Unfortunately, this cannot be 
done in any meaningful way for the atomic and ionized me- 
dia, whose true densities are not known. In contrast, crude 
upper limits can be derived for the two components of the 
molecular medium. Remember that the filling factor is by 
definition the ratio of space-averaged to true density. The 
H2 space-averaged density attains its maximum value in the 
CMZ, where (nH2)max = 150 cm~^ (see Equation [T5|) . If, 
as argued by Rodriguez-Fernandez et al. (2001), ~ 30% of 
the interstellar molecular gas resides in the warm compo- 
nent (see section [221) ^nd if the cold and warm components 
have true densities ~ lO'^'^"'' cm~^ and ~ 10^'^ cm~'^, re- 
spectively, then their filling factors are < 3.3% and < 14%, 
respectively. 
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